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ABSTRACT 


The detailed design of a beam-powered transatmospheric vehicle, “The Apollo Lightcraft ”, 
was selected as the project for the design course. The vehicle has a lift-off gross weight of about six 
(6) metric tons and has the capability to transport 500 kg of payload (five people plus spacesuits) 
to low Earth orbit. Beam power was limited to 10 gigawatts. The principle goal of this project is to 
reduce the LEO payload delivery cost by at least three orders of magnitude below the Space Shuttle 
Orbiter - in the post 2020 era. 

The completely reusable, single-stage-to-orbit shuttlecraft will take off and land vertically, 
and have a reentry heat shield integrated with its lower surface - much like the Apollo Command 
Module. At appropriate points along the launch trajectory, the combined-cycle propulsion system 
will transition through three or four airbreathing modes, and finally use a pure rocket mode for 
orbital insertion. 

As with any revolutionary flight vehicle, engine development must proceed first. Hence, the 
objective for the Spring semester propulsion course was to design and perform a detailed theoretical 
analysis on an advanced combined-cycle engine suitable for the Apollo Lightcraft. The class deter- 
mined that only three airbreathing cycles would satisfy the mission, and that the ramjet cycle is 
unnecessary. 

The preliminary theoretical analysis of this combined- cycle engine is now complete, and 
the acceleration performance along representative orbital trajectories has been simulated. Average 
vehicle acceleration is approximately 4-5 G‘s. Transition between engine modes occurs at Mach 3, 
11 and 25+. Beam power can be reduced to as low as 2.5 billion Watts without sacrificing vehicle 
performance. The Lffi propellant requirement is typically 300 kg, or roughly 5. % of the vehicle 
lift-off weight, for delivery to a 100 nautical mile orbit. The total beam energy requirement is 520. 
GW-sec for this boost mission: 

Therefore, the propellant cost at current bulk LHi rates is $ 975: the energy cost is 1 2455. 
assuming present wholesale hydroelectric power rates. The total round trip cost is only $ 3430, or $ 
686 per person. This represents a payload delivery cost of $ 3.11/lb, which is a factor of 1000 below 
the STS. 

The Apollo Lightcraft concept is now ready for a more detailed investigation during the 
Fall semester “Transatmospheric Vehicle Design "course. The class will divide itself into smaller 
design groups (5 students in each) dedicated to: a) Aerodynamics, b) Propulsion, c) Structures, d) 
Thermal Analysis, e) Flight Control Systems, f ) Optimal Trajectory Analysis, g) Human Factors and 
Life Support Systems, and h) Power-beaming Architecture These smaller design groups will interact 
with each other, the TA (a graduate student), and the instructor to finally evolve an integrated 
conceptual design for the Apollo Lightcraft vehicle. 
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CHAPTER I 
INTRODUCTION 


Contained herein is the preliminary design and performance analysis of an advanced com- 
bined cycle engine suitable for application in a small reusable single-stage-to-orbit (SSTO) shuttle- 
craft. The innovative combined-cycle engine is designed to utilize highly-energetic beamed-energy 
sources (laser, microwave), and propell a manned transatmospheric vehicle based on the Apollo 
Command Module (CM). 

Dubbed the “Apollo Lightcraft ", this vehicle must transport a five-person crew (total 
payload = 500 Kg) into low Earth orbit (LEO) in three minutes, or anywhere on the globe in 
one half hour. Fully loaded for launch, the Lightcraft would weigh about six metric tons which, 
incidently, is equal to the dry weight of a Harrier jumpjet, Learjet, or thirty-foot long Winnebago 
motor home. 

This study attempts to “create the blue-prints "for revolutionary globe-trotting, laser-riding 
shuttlecraft of the 21st century. The design group hopes that such visionary approaches will finally 
enable large scale access to space by reducing the payload transport costs by a factor of 1000 below 
the Space Shuttle Orbiter. It is with this bold objective in mind, that the present team set out to 
prove the feasibility of the “Apollo Lightcraft ". 

A. Vehicle/Propulsion System Description 

Figure 1 is a sketch of a 10 GW orbital shuttlecraft that would use advanced combined-cycle 
engines for propulsion. The man/machine interface is inspired by the Apollo space-capsule (see Fig. 
2). It is interesting to note that this minimum weight/volume spacecraft approach was originally 
taken at a time when chemical rocket propulsion systems were in their infancy and rapidly maturing 
toward larger launch-mass (payload) capability. 

As indicated in Fig. 3, the Apollo Lightcraft would have a diameter of 5 m, gross weight 
at takeoff of 5550 Kg and payload of 500 Kg. The five person crew would be propelled to orbit in a 
single-stage vehicle that has a reentry heatshield integrated with its entire afterbody. As portrayed 
in Fig. 4, the annular shroud would be translated to the full forward position, and movable thermal 
protection system (TPS) tiles would slide forward to close off the annular gap (i.e., between the 
vehicle centerbody and the shroud), just prior to reentry. 

Pictured in Fig. 5 is one concept for a 10 GW Satellite Solar Power Station (SSPS), an idea 
first conceived by Peter Glaser of Arthur D. Little, Inc. (Cambridge, Ma.). Now imagine a hundred 
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of these SSPS’s linked together into a space-based “power grid "in which “wireless "transmission of 
energy is accomplished by laser (or microwave ) beams - sometime in the early 21st century. The 
succesful mating of advanced beam-powered flight vehicles (such as the Apollo Lightcraft) to highly 
reliable lasers and satellite solar power stations would revolutionize air and space transport as we 
know it today. Clearly, all necessary safeguards must be designed into the space power system - 
in order to provide extreme levels of network security, and accountability for every Joule of energy 
beamed throughout cislunar space. 

As indicated in Fig. 3, beam power is first received by the Apollo Lightcraft across a 4.25 
meter diameter centerbody primary optics at an intensity below that which would ignite a high 
temperature plasma (e.g., 1 - 2xl0 4 W fern 2 for 10.6 micron light). This plasma-ignition threshold 
is lower for longer wavelengths, and higher for shorter wavelengths. Hence, the greatest power is 
delivered at the ultraviolet limit of about 0.35 pm - beyond which atmospheric scattering dominates. 
The primary optic surface area of 8 x 10 4 cm 3 (see Fig. 3) would permit the reception of 0.8 to 1.6 
GW of continuous laser power at 10.6 pm. At the ultraviolet limit, this value can be increased by 
a factor of five or so. It should be noted, however, that the first two engine modes for the Apollo 
Lightcraft are pulsed, and one must be careful to reduce this peak continuous power level by the 
effective duty cycle. 

The Apollo Lightcraft combined-cycle engine would permit vertical takeoffs and (powered) 
landings. While accelerating to orbit this engine would “shift gears ”, much like a 4 or 5 speed high 
performance sports car. The advanced engine would have three or four airbreathing modes, plus a 
rocket mode used for orbit insertion. These engine cycles (or modes) are chosen to maximize flight 
performance throughout a given range of Mach numbers, with distinctly identified transition points. 
It is evident that the choice of specific combined-cycle engine schemes has a large effect upon the 
overall vehicle configuration. 

For the Apollo Lightcraft, the most logical choice for the combined- cycle engine is as follows: 

Mode 1 - ERH thruster; 

Mode 2 - ramjet (optional); 

Mode 3 - scramjet; 

Mode 4 - MHD - fanjet; and, 

Mode 5 - rocket. All engine modes are powered by beamed energy. 

Fig. 5 portrays the Apollo Lightcraft in the External-Radiation- Heated (ERH) thruster 
mode, which is used for powered vertical takeoffs and landings, as well as acceleration runs up to 
Mach 3. Note that the ERH thruster reaction surface (see Fig. 6) is easily combined with a reentry 
heat shield - with little additional weight penalty. 
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It has been discovered that a “pop-up " maneuver (using the ERH thruster) must be utilized 
before starting the acceleration run to orbit. Under a low power microwave beam link, the shut- 
tlecraft could climb vertically through any weather (clouds, fog, rain, snow, etc.), until an altitude 
of 20 - 40 kft. is reached. Then the Lightcraft would rotate to engage high power laser beam link 
from a low altitude relay station (e.g., 185 km), and start accelerating. Not only does this procedure 
guarantee maximum acceleration performance, but it also insures all-weather capability. 

Upon reaching Mach 3, propulsive power would be transferred to the ramjet (or perhaps 
directly to the scramjet) mode. The incoming laser power would be absorbed within the annular 
shroud region. An annular “planar heater ” (pictured in Fig. 7) would be created by repetitively 
pulsed LSD waves that propagate at right angles across the duct flow. Note that the planar heater 
would promote turbulent mixing of the internal flow and enable an exceptionally short shroud. 

As the vehicle continues to accelerate, this duct flow would increase to supersonic velocity, 
and the engine would automatically transition to the scramjet mode. As indicated in Fig. 8, the 
annular shroud would translate aft with increasing flight Mach number such that the conical bow 
shock is always attached at the shroud forward lip (i.e., beyond Mach 3, the design Mach number 
of the inlet). 

At some point beyond Mach 11, frozen flow losses would dominate, and net thrust would 
fall to zero. The Apollo Lightcraft engine would then transition to the MHD-fanjet propulsion cycle 
portrayed in Fig. 9. As indicated, laser-heated H 2 rocket-gas generators would extract 40 to 50 % 
of the enthalpy as electric power and deliver it to an airbreathing “electric fan " (see Fig. 10). 

The Apollo Lightcraft annular duct would be converted to an MHD air accelerator (see Fig. 
11) simply by energizing a series of cryogenic or superconducting electromagnets . As shown in Fig. 
12, these electromagnets are integrated within the shroud support struts and generate a confined 
toroidal magnetic field. The magnetic field configuration is designed to close upon itself, such that 
the field lines do not arc out into the surrounding airspace external to the engine. It is important to 
note that the MHD fanjet can exhibit coupling coefficients in the same range as a pure laser-heated 
rocket (i.e., 50 - 100 N/MW), but can have specific impulses an order of magnitude higher (e.g., 
10,000 to 20,000 sec., vs. 1000 to 2000 sec.). 

Once orbital velocity has been reached at an altitude between 200 to 250 kft., the MHD- 
fanjet cycle would be transitioned to pure rocket mode (i.e., electromagnets are shut off), and the 
vehicle would pitch up to leave the atmosphere. The twelve short laser heated rocket exhaust nozzles 
are positioned radially about the vehicle afterbody - which acts as a plug nozzle - to secure a large 
expansion ratio at high altitudes. The mass flow rate of hydrogen would roughly double to ll.Kg/s, 
and the stagnation temperature would fall to perhaps 8000 K, with the stagnation pressure reduced 
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Figure 9. Diagram of MHD-fanjet cycle (an electric air-turborocket). 
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Figure 10. Airbreathing MHD accelerator concept (hypersonic mode). 
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Figure 12. Hypersonic MHD accelerator concept (view from rear, looking forward). 
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to 60 bare. Under these conditions, the rockets would develop a specific impulse of 2000 to 2400 
s. With a thermal efficiency of 66 %, the engine could produce a coupling coefficient of 67 N/MW. 
At a vehicle mass of 5560 Kg, this gross thrust level of 2.35 x 10 6 N would result in a maximum 
acceleration of 4.3 g (i.e., if vehicle drag is neglected). Since the vehicle would have depleted some 
fraction (e.g., one-third), of the onboard propellant by this time, the acceleration could be a little 
higher. If a lower acceleration is desired, it could be exchanged for a higher I,p from these rocket 
engines. 

Fig. 13 attempts to summarize in one chart the coupling coefficients anticipated from laser- 
heated rockets, in comparison with most known airbreathing engine cycles. The goal of advanced 
trans-atmospheric “combined-cycle ” engines is, of course, to reduce the fuel fraction required to 
carry a given payload into orbit. Indications are that laser- powered, SSTO combined-cycle engines 
could enable fuel fractions as low as 5 to 10 %. Saying it another way, mass fractions (initial/final) 
of 1.11 or less may be altogether feasible! 

Once the decision is made to leave low Earth orbit, the Apollo Lightcraft rocket engines 
would be fired (briefly) to slow the vehicle by perhaps 100 m/sec. The spacecraft would then reenter 
the atmosphere, shielded by a thermal protection system similar in design to that used on the earlier 
Apollo capsule. Eventually the Lightcraft would decelerate to a specific subsonic terminal velocity 
(i.e., a free-fall), and just prior to landing, the ERH thruster would be fired for a few seconds to 
provide a braking force. The long-term goal for such advanced laser-powered shuttlecraft should be 
the same kind of operational flexibility demonstrated by today’s airlines, or better. 

Since Lightcraft will have VTOL capability, they could in principle be set down almost 
anywhere. They do not require a 10,000 foot long runway; just a helicopter pad will do. Also, 
since the power beam provides the energy for propulsion, the LH 2 propellant (used during the high 
performance MHD-Fanjet and rocket modes) might possibly be replaced by something as simple 
as water! Although the acceleration performance and specific impulse would suffer, the cost of de- 
ionized water is negligible relative to LH?. Hence the overall payload delivery cost might be roughly 
the same - and you could “fill’er up " at your house with the garden hose! Finally, judging on the 
basis of the available tankage volume in the Apollo Lightcraft, the total “fuel load " might account 
for 6 to 46 % of the launch weight, depending on the performance of the combined-cycle engine 
scheme and the choice of propellant (ie., LHi vs. H 2 O, respectively). 

In some post 2020 era, space travel will be commonplace, and the present large “standing 
army * ground crew will be replaced by efficient macro-computers - which schedule launch windows, 
effortlessly deliver beam power from the Space Power Grid to hundreds of spacecraft simultaneously, 
and electronically mail out the end-of-the-month billings to millions of satisfied “jet-setters ". The 
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dominant costs will be for propellant and power. Payload delivery costs will eventually plunge a 
factor of 1000 below that of the Space Shuttle Orbiter. 

B. Design Study and Final Report Orgnalsatlon 

The Apollo Lightcraft combined-cycle engine study was carried out as a “proof-of-concept 
” investigation, and consisted of the following eight tasks: 

Task 1 : Compute Maximum Available Beam Power vs. Wavelength 

Task 2 : Select Engine Modes and Transition Points 

Task 3 : Develop Detailed Analytical Models for Engine Performance 

Task 4 : Conceive Integrated Engine/Vehicle/Optics Configuration 

Task 5 : Determine Forebody, Shroud and Base Drag vs. Mach No. and Altitude 

Task 6 : Analyze Altitude vs. Mach Number Performance of all Engine Modes 

Task 7 : Estimate Shuttlecraft Mass Breakdown 

Task 8 : Predict Shuttlecraft Performance Capability Along Orbital Trajectory 
Principal results of the first six tasks are reported in Chapts. II through V of the present 
report. Chapter II presents the model for the ERH thruster, chapter III, the Ramjet/Scramjet 
model; chapter IV, the MHD-Fanjet model, and finally chapter V, the rocket model. 

These detailed analytical models were developed for assessing the thermodynamics and 
gasdynamics of all four airbreathing engine modes given specific data on component technology. 
The models were assembled to facilitate the computation of thrust, coupling coefficient, specific 
impulse (7*p), and beam power - as a function of flight conditions (i.e., Mach number and altitude). 
These models permit the computation of gas stagnation and static conditions, cross-sectional flow 
dimensions and flow Mach number at all important engine stations. Complete engine diagrams 
showing the important engine stations are presented with the performance in later chapters of this 
report. 

The importance of inlet centerbody cone angle, shroud configuration and plug nozzle design 
upon engine total pressure recovery were found to be significant. The choice of inlet cone angle, 
for example, is a compromise between the needs of two separate engine modes: the MHD-fanjet 
benefits from a blunt cone inlet to enhance flow ionization; whereas, the ramjet and scramjet would 
like sharp cone inlets for maximum total pressure recovery. Note also that light-weight structural 
design dictates minimal enclosed volume, which points to the blunt cone solution. 

Using the detailed analytical models developed in Task 3 and the data on specific engine 
components, the altitude/velocity performance was predicted for the airbreathing comined-cycle 
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engine. Partial (throttled) and full thruat in each mode were characterized. Optimum performance 
in all modes was determined, as well as the most favorable transition Mach number between modes. 
Detailed computer models were used to carry out these calculations. Principal results included 
graphical plots of thrust, specific impulse, coupling coefficient, propulsive efficiency and beam power 
vs. flight Mach number - with altitude as the parametric variable. 

The results of Task 8, which predicted the performance capability of the Apollo Lightcraft, 
are presented in chapter VI. In addition to determining the acceleration schedule (i.e., G-level vs. 
time) and the “time-to-climb ” to orbit, a primary output of this task was to predict the overall 
vehicle mass ratio (i.e., final-to-initial weight), time-average beam power along the trajectory, and 
total boost energy in gigawatt-seconds required to deliver the 500 Kg payload to orbit. The group 
assumed a beamed energy cost of 1.7 cents/KW-hr (current NE Canadian wholesale hydro-power 
rate), and a price Of $ 3.252/Kg for LHi (current Shuttle Orbiter rate). 

Finally, chapter VII summarizes the principal conclusions reached in this initial proof-of- 
concept study on the Apollo Lightcraft combined-cycle engine. 
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CHAPTER H 
ERH THRUSTER 


The External Radiation Heated (ERH) thruster, as its name implies, is an external air- 
breathing engine which utilizes a high intensity laser beam to produce thrust. The thruster itself is 
extremely simple, consisting only of an impulse surface, which is integrated in the aft section of a 
flight platform and able to withstand high temperatures and pressures. In the present analysis, the 
ERH thruster surface is a circular frustum which makes up the aft section of the Apollo Lightcraft, 
see Figure 14. To generate thrust, a series of laser-induced cylindrical blast waves are initiated 
adjacent to the thruster surface; as these blast waves expand, impulse is delivered to this surface. 
To understand the principles of this thruster, it is first instructive to examine the laser-induced blast 
wave phenomenon which is fundamental to its operation. 

A. BASIC OPERATING PRINCIPLES 

When a high intensity laser beam is focused, a detonation wave forms: the structure of 
which consists of a shock followed by an absorption zone. Within the absorption zone, laser energy is 
transferred into the fluid, in the form of thermal energy, by the mechanism of inverse bremsstrahlung. 
This Laser-Supported Detonation (LSD) wave propagates up the laser beam at supersonic velocities 
and leaves behind a high energy, high pressure plasma zone. This, in turn, expands into the ambient 
air much like a blast wave generated by a conventional chemical explosion. Because LSD waves 
propagate toward the laser source much faster than the radial expansion of the hot plasma behind 
the detonation front, the resulting blast wave is cylindrical in shape. 

The ERH thruster utilizing LSD blast waves has the following operation. First, several 
laser beams projected parallel to and across the thruster surface, are focused to initiate LSD waves. 
As these LSD waves propagate across the thruster plate, high pressure cylindrical blast waves are 
formed and subsequently expand. The portion of the blast wave in contact with the thruster surface 
exerts a pressure force upon the surface. It is this impulse loading that results in thrust. Once the 
Blast wave pressure decays to the local ambient pressure, the thruster surface must be cleared of the 
hot plasma and replaced with ambient air before another set of LSD waves can be initiated. This 
exchange of the hot plasma with cool, unprocessed air is equivalent to the heat rejection portion of 
a conventional thermodynamic cycle. In the present analysis, this exchange process is referred to as 
surface refresh. 

From the above discussion, it is evident that the ERH thruster is a repetitively-pulsed 
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engine, in which its Pulse Repetition Frequency (PRF) is governed essentially by the blast wave 
decay time and the thruster surface refresh time; the total cycle time is dependent on the local 
ambient air pressure and density. This implies that thruster performance varies with altitude and 
flight velocity. Before the ERH thruster performance is analyzed, the structure of the LSD wave is 
investigated and a static model for the ERH thruster is developed. 

B. LSD BLAST WAVE ANALYSIS 

A laser supported detonation (LSD) wave can be analyzed as if it were a chemical detonation; 
the chemical heat of reaction is replaced by the absorbed laser heat flux. Therefore, the ZND 
detonation wavemodel, developed by Zel’dovich, Von Newman, and Doringfor chemical detonations, 
can be adapted and used to analyze the LSD wave structure. In the ZND model, chemical detonation 
wave structure is modeled as a one-dimensional shock front followed by a high speed deflagation; for 
LSD waves the deflagration is replaced by a thin absorption zone. In LSD waves there is a coupling 
between the fluid dynamics and the absorption kinetics. The leading shock ionizes the air thereby 
allowing the formation of the absorption zone, and likewise, the laser energy absorbed within this 
zone drives the shock. Because of this coupling, the detonation wave structure occurs only when 
laser intensities are above 10® - 10 7 W/CM i (for 10.6 pm radiation). Raizer(13], using the ZND 
wave structure, was able to derive the following wave velocity expression for a LSD wave propagating 
into quiescent gas. 


Vo - (2(t> - 1)±!‘' S (II -1) 

P o 

Up to this point, the detonation wave has been viewed from an Eulerian frame of reference It 
is, however, generally easier to examine the detonation wave structure from a Lagrangian frame. In 
the Lagrangian frame, the observer is assumed to be traveling on the shock front, typically referred 
to as shock coordinates. To convert from Eulerian coordinates to shock coordinates the following 
transformation is applied: 

« = «' - Vo (//- 2) 

In Fig. 15 the structure of the LSD wave is shown in shock coordinates. Observe that a radial 
rarefaction fan begins to propagate inward toward the blast wave center directly behind the absorp- 
tion zone. Also note, that regions of steady pressure and expansion velocity exist ahead and behind 
the rarefaction fan. In an actual blast wave, the converging expansion fan is reflected at wave’s 
centerline and forms a second rarefaction wave which propagates outward from the wave’s center. 
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In this analysis, only the first rarefaction wave will bo of interest. Once the first fan reaches 
the centerline, at location X' in Figure 15, the blast wave decay is approximated using self-similar 
theory. Previous investigators (7, 12] have sucessfully applied self-similar blast wave theory away 
from the LSD detonation front. Here, the blast wave will be analyzed in two sections: first, using 
the method of characteristics in the region between the absorption front and the plane where the 
first fan reaches the center line; and thereafter, using self-similar blast wave theory. 

In this model, axial relaxation of the LSD blast wave is neglected, and radial expansion is 
assumed to dominate the pressure decay process. Also, the axial velocity is assumed small when 
viewed from the Eulerian reference frame. These assumptions permit the use of the following one- 
dimensional unsteady flow equations to describe the blast wave expansion: 


, dp . du , pu rt 

W + , W + 'W + r 0 

{II - 3a) 


1 

co 

o- 


{I I -3c) 


For the present development, the last term of the continuity equation is dropped. Using an order of 
magnitude analysis, it can be shown, that the error introduced by neglecting this terms is small; its 
removal simplifies the development of the characteristic equations. Using the method of characteris- 
tics, the above system of partial differential equations are recast into ordinary differential equations, 
such that the fluid properties are described along some characteristic path. Hence, the continuity 
and momentum equations can be combined and manipulated to give the following result: 

« + -~r = «/ + = Const. {II ~ 4o) 

7-1 

« - = J~= Const. (II - 46) 

7-1 

These equations are the well known Riemann invariants which are valid along the following paths 
in the length-time plane. 

— - a it + C ; C + Characteristic (// - 5a) 

at 

^ s u - C ; C~ Characteristic {II - 56) 

at 

The above four expressions are fundamental to the analysis of LSD blast waves near the detonation 
front. 

As previously mentioned, the pressure decay process for the rest of the blast wave is ap- 
proximated using seif-similar cylindrical blast wave theory. Expressions predicting this type of wave 
decay are easily derived from dimensional analysis. 
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Assume that the column of high pressure plasma generated by a LSD wave can be approx* 
imated by a cylindrical blast wave resulting from an intense line source explosion. See Figure 16. 
The dependent variables of this blast wave expansion are velocity (U), pressure (P), and density (p), 
while the independent variables are blast wave radius (R), and time (t). The physical constants of 
this decay problem are the energy released during the explosion, as well as, the initial density and 
pressure. The dimensions of these variables are: 


{®}=™ 

(II -6a) 

M = M 

(I I -6b) 

II 

(I I -6c) 


(II- 6d) 

It can be shown, that the initial pressure is not important in strong shock waves; thus, the only 
parameter constants are the energy released and the initial air density. The objective of dimensional 
analysis is to form nondimensional parameter groups. Observe that the two independent variables 
can be combined to form the following similarity parameter 

"° = (Z) 1/5*17* = Const • 

(I I -7) 


Hence, r is proportional to t -1 / 2 . Using the strong shock assumption, following expression can be 
obtained: 

r *^)*’* 2 («-«) 

Differentiating Eqn. 7 and substituting it into the above expression results in an expression for 
pressure: 

P = (”-9) 

In the above equation the bracketed terms are constant and it is possible to form the following ratio: 

j£rO'‘ 

An additional ratio, can be derived, using Eqn. 7: 



( C 7 )J 


(JJ-U) 


Hence, if blast wave conditions are known for a given time, then the above expressions determine 
subsequent blast wave pressure decay. 
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To determine the fluid conditions directly behind the absorption zone, in the Chapman- 
Jouguet plane, the following equations developed for the ZND model can be used: 


« CJ = 


_ Vd 
7+1 


(// - 12 a) 




(//- 126) 

C ? c/® 7 «cj (//- 12 c) 

Recall that the Chapman-Jouguet plane is located just behind the reaction zone; in the case of LSD 
waves it is behind the absorption zone. Chapman-Jouguet theory predicts that stable detonation 
waves exist only when the relative velocity of the flow in the Chapman-Jouguet plane, with respect 
to the wavefront, is equal to the local sound speed; therefore: 

V D -ucj = C C j = \^-]i (//- 13) 


CJ 


Combining the definition of the ideal gas sound speed, the expression for polytropic expansion and 
Eqns. 12c and 13 the pressure behind the detonation is given by: 




{II - 14) 

7 + 1 7 + 1 ,v v ’ 

Assuming that the axial velocity is zero, the initial impulse pressure and sound speed are given by 
the following expressions, respectively: 




(7/ -15) 


r - v d 
Cl ~~2 


{II - 16 ) 

Using the above results, the Riemann invariants and the self- similar expansion expressions, it is 
now possible to determine the expansion of the blast wave and subsequent thrust generation. The 
pressure decay near the detonation front is now examined. 

From Fig. 17 it is seen that the first rarefaction fan is a simple backward facing wave; thus, 
flow conditions behind it can be calculated using the invarient along the C+ characteristic. Because 
the radial velocity ahead of the expansion fan is equal to zero, the invariant is written as: 


2 _ 2 ^ 

__C l = t fc + _C 2 

Combining this with the perfect gas relations, the above equation can be written as: 


{II ~ 17) 




{II - 18a) 
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(II-lSb) 


= -(£>»(£)*! 


("-«•) 

This expression can then be solved numerically with the following well established shock tube result, 
to determine the pressure behind the rarefaction fan (Pm). 


<» “<!?>♦, I -‘,1 


2 . 

ho + 1)$*+ (7o - 1) 


I* 


(// - 19) 


Substituting this result into equation 18c, the particle velocity behind the refraction fan can be 
determeined. Knowing the particle velocity, it is now possible to find the location of the radial 
shock front. At X\ the radius of the blast wave is defined as Rr«j. This reference radius is used to 
start the self-similar analysis. 


C. STATIC THRUST MODEL 


The impulse loading on the thruster surface due to an expanding cylindrical plasma blast 
wave, is given by the integral relation: 


I 



PiAdt 


(II - 20 ) 


Using the above expression, the impulse delivered to the thruster plate can now be found. 
To determine the impulse contribution of Region 1 (see Fig 15) note that its geometry 


remains fixed as the LSD wave propagates over the thruster surface. The location of X' can be 


determined by integrating Eqn. 5b, solving for the integration constant and multiplying the result 


byVb. 


X' = 2 R, 

Therefore the impulse contribution of this region is simply: 


(/I -21) 


/i=2r.(r f Pi+r re /P 2 )f (11-22) 

where the t is equal to the time necessary for the detonation front to propagate across the 
thruster surface. Since the laser pulse is terminated when the detonation front reaches the end of 
the thruster plate, t is equal to the laser pulse duration time t p . 
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The impulse contribution for the region X > X' in Fig. 18, can be determined from 
equations 10, 11 and 20, which are combined to give: 

h = r rt ,PrefAXi f° = 2 Ifre/Pre/WK^ 2 -)* - IjA* (// - 23) 

Jt r ,f W ‘re/ 

where to is the time necessary for the blast wave to decay to the local ambient pressure, and t re { is 
the time required for the first rarefaction fan to reach the blast wave center. To account for varying 
flow properties along the plug nozzle, the impulse in the radial expansions regime is numerically 
integrated over the thruster surface, using a given step size of delta x. Thus, the total impulse 
generated by the LSD wave is the sum of the contributions from the two expansion regimes, as given 
below: 

/total = /i + ^ /< (II- 24) 

•si 

Note that the impulse force acts along the thruster’s surface normal. For the ERH thruster engine 
configuration examined, all of the impulse loading predicted by the above expression does not con- 
tribute to engine thrust; only the impulse vector component aligned with the vehicle flight direction. 
Thus, magnitude of the impulse component resulting in upward thrust is: 

/usable = /total cos a (II- 25) 


where a is the base angle of the truncated conical plug. 

In trajectory calculations is convenient to use time averaged thrust. This mean thrust is 
defined as: 

f Woven _ /TOTAL 


7 = 


1 


tcYCLE 


1 


( 11 - 26 ) 


t CYCLE 

where I is the impulse delivered per laser pulse, and T c j C u is the blast wave period. The blast wave 
period is simply the total expansion time of the LSD blast wave, plus the refresh time; as given by: 


f CYCLE = 


1 

TW 


= to + ^REFRESH 


( 11 - 27 ) 


The refresh time is the time necessary for the hot expanded gas to clear the thruster surface and 
be replaced by the ambient air. Currently, an adequate model for the refresh process does not exist 
and in this analysis the refresh time is set equal to zero. 

Input energy per laser pulse is evaluated as follows: 



(II - 28) 


In the derivation of this equation it is assumed that the laser beam has a semicircular shape and 
uniform intensity. Recall that the laser energy is absorbed in a thin region behind the detonation 
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front. Thus, t p is the time requried for this LSD wave to transverse the ERH thruster plate. To 
convert the pulse energy to time average laser power the following equation is used: 

^AVB = -r—T (It - 29 ) 

►CYCLE 

Finally, an important measure of thruster efficiency is the coupling coefficient, which is the ratio of 
usable single- pulse impulse to laser input energy per pulse. Therefore, the coupling coefficient is 


defined by: 


nn _ Aj SABLE _ 

cc I? 



(7/ -30) 


D. DYNAMIC THRUSTER MODEL 

The above analytical model does not consider the effects of flight platform velocity on ERH 
thruster performance. This model is sufficient for static and low subsonic flight velocities, but when 
flight velocities are supersonic, ram drag, compressiblity effects and flow expansion over the aft plug 
nozzle (the thruster surface) become important. For an airbreathing engine net thrust is the force 
resulting from the addition of all pressure and viscous forces, excluding the external drag forces; 
that is, gross thrust minus ram drag. Note external drag forces caused by the "installed” engine are 
typically included with vehicle drag. As a result, it was necessary to develop a dynamic model for 
ERH thruster performance. 

In supersonic flight, an attached oblique shock forms at the conical spike tip. The air is 
compressed after passing through this oblique shock, and it is turned to flow tangent to the forebody 
surface; the entire forebody of the flight platform is essentially an isentropic spike inlet. The annular 
cowl, which circumscribes the craft’s midsection, then redirects the flow parallel to the thruster 
surface. In the ERH thruster operation, the cowl only turns the flow, and does not decelerate or 
accelerate it. Additionally, the cowl is assumed to have a sufficiently large capture area such that no 
spillage occurs beyond the design mach number of three. With these cowl characteristics, no normal 
shock forms at the cowl entrance and the exit mach number is equal to that at the cowl entrance. 
Upon leaving the cowl, the ducted air expands along the thruster surface as if it where emerging 
from a free expansion type plug nozzle. 

Therefore the only stagnation significant pressure loss occurs as the air passes through the 
oblique shock. To calculate such losses, an experimentally determined pressure recovery inlet sched- 
ule developed by Marquart|l5] for the XRJ59-MA-3 inlet was used. In the present analysis, it is also 
assumed that the air flow through the cowl remains undisturbed by the LSD wave propagation and 
subsequent blast wave expansion. Such an assumption appears to be justified, when the maximum 
expansion of the blast wave is small compared to cross sectional flow area of the cowl. 
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In supersonic flight, pressure variations along the thruster surface become pronounced and 
must be considered when calculating the LSD blast wave expansion. When calculating this pressure 
variation, it is assumed that the compressed air flow leaving the cowl expands as shown in Figure 
19. To calculate the pressure for a given thruster surface location (s) it is necessary to determine 
the local cross sectional area of the expanding gases. 

Typically, for a non-truncated plug nozzle, the pressure at the end of the plug is always 
equal to Po, the ambient pressure. Using this condition and assuming isentropic flow over the plug 
the following expression can be solved for the mach number at the plug end (A/ e ). 

= (7/ -31) 

where Pu« denotes the stagnation pressure at the cowl exit. As previously mentioned this is deter- 
mine using the Marquardt inlet data. Once M e is known, the cross-sectional flow area, A e , can be 
obtained using: 

Aee _ jWsi 

At Met 

where A«« and Me* are the cross-sectional flow area and the mach number at the cowl exit respec- 
tively. 

From the nozzle’s geometry: 

R* = \J^ (// - 33) 

If the cross-sectional flow area is assumed to vary linearly along the thruster surface, then: 



A * t\R% - (— )I 
1 * v cosa ,J 


(J/-34) 


In this equation, a is the plug base angle and s denotes an arbitrary position along the plug surface. 
For a given s, the coresponding mach number can be found implicitly using: 


Now the local expansion pressure can be calculated from: 


P(>) = 




(II - 35) 


(II - 36) 


In the dynamic model, the LSD blast waves must expand to this pressure. 

To calculate the net thrust, a control volume is constructed as shown in Fig. 20. By choosing 
the control volume boundary adjacent to ERH thruster surface, the complex flow conditions due 
to the interaction between the blast wave expansion and the air flow directed by the cowl can be 
neglected; this simplifies the drag analysis. 
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Figure 



Recall that in the development of the LSD blast wave model, the ambient pressure was 
already subtracted from the integrated impulse, hence it is necessary only to consider the pressure 
difference across the cowl. The verticle force due to this pressure difference is given by: 

Pressure Drag = (Po - Pee)A c « sin a (II - 37) 

Note that there is a force due to the momentum change of the incoming air flow is redirected to 
convect over the thruster surface. This force is calculated using the following expression: 

Cowl Drag = sin a (II- 38) 

Therefore the net thrust is given by: 

Tnct = T- Pressure Drag - Cowl Drag (II - 39) 


E. RESULTS 

Performance maps for the ERH thruster engine were generated by incorporating the pre- 
ceding static and dynamic models into a comprehensive computer simulation. Using this computer 
model, the effects of flight mach number and altitude on net engine thrust, as well as several other 
performance parameters was investigated. 

For the Apollo Lightcraft, the following ERH engine configuration was assumed. As previ- 
ously mentioned, the thruster surface is a truncated conical plug which also comprises the entire aft 
section of the vehicle. This plug has a semi-vertex angle of 45 degrees. During ERH thruster opera- 
tion, 48 LSD waves are ignited; they are equally spaced around the circumference of the plug. The 
laser beams which initiate these waves are assumed to provide a uniform intensity of 5rl0 8 W/cm 2 , 
over the semi-circular spot of 0.5 cm radius. 

In Fig. 21 net thrust variation is presented as a function of flight mach number and altitude. 
Recall that net thrust is the time-averaged thrust minus ram drag. For a given altitude, observe 
that thrust drops off at some critical mach number. This is due, in part, to the increased ram drag 
experienced by the ERH thruster engine at higher mach numbers. 

Fig. 21 also shows that mean thrust decreases with increasing altitude. This result is not 
unexpected since at higher altitudes atmospheric pressure and density decrease. The ERH thruster, 
like other airbreathing engines, relies on momentum transfer of incoming ambient air to generate 
thrust; therefore, as the density decreases, so does engine thrust. 
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Figure 21. ERH thruster performance - Time average net thrust. Figure 22. ERH thruster performance - PRF. 


Another factor which causes a reduction in net thrust, for this ERH engine configuration, is 
illustrated in Fig 22: the maximum pulse repetition frequency (PRF) rate. Note that PRF decreases 
with increasing altitude. During the ERH thruster operation, LSD blast wave pressure must expand 
to local ambient pressure before the thruster surface can be refreshed and new LSD waves initiated. 
At higher altitudes, the local static pressure is lower, so the blast wave expansion time (to) is longer. 
Since the high-internal pressure of the plasma blast wave decreases rapidly as the wave expands, 
most of the engine thrust is generated during the initial blast wave decay. Therefore, an increase 
in the LSD blast wave expansion time (to) does not greatly increase the impulse delivery to the 
thruster plate, and time averaged thrust must decrease. 

Closely related to the PRF is the time-averaged laser power, its variation with altitude and 
flight mach number can be seen in Fig. 23. Also, in Fig. 24, the predicted coupling coefficient 
performance is plotted, with the dotted line indicating the theoretical limit for coupling coefficient 
of airforeathing engines. Since the calculated coefficients at high altitudes (e.g., 30km) so closely 
approach the theoretical limit, it is suspected that this first order analysis of ERH engine performance 
may be overly optimistic. There are a number of reasons for these high values. 

In the current model, as mentioned above, the refresh time (t r «fr«ah) is set to zero. In an 
actual engine, t r «frMh would have some finite value, which would cause a decrease in the PRF and 
hence, the time averaged thrust. A second model simplification, which probably contributes to 
producing optimistic results, is the omission of axial rarefaction waves within the LSD blast wave 
simulation. Inclusion of axial expansion effects in this blast wave model, would result in a reduction 
in the impulse delivery to the thruster surface. At the present time, the combined effects of these 
simplifications on ERH thruster performance have yet to be quantified. Therefore, to be conservative 
the net thrust values used in the trajectory analysis are reduced by 33%. 
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CHAPTER m 


RAMJET/SCRAMJET MODE 

The scramjet engine will begin operation at about Mach 3.0 and 20000ft, where the ERH 
thruster cuts off. The scramjet engine produces thrust by adding heat to supersonic internal engine 
dr. In current scramjets, liquid chemical fuel is burned to heat the engine working fluid. In the 
Lightcraft engine, however, energy is added by focusing a laser beam into the combustion chamber 
(see Fig. 25). The flow is heated to high temperatures by passing through a “planar heater " 
produced by a repetitively-pulsed LSD wave which propagates radially outward from the centerbody 
to the cowl. The laser is shut off when the LSD wave reaches the cowl inner lip, where the secondary 
optics are located. 

A. INLET SPIKE 

Since the scramjet is used in supersonic/hypersonic flight from about Mach 3.0 to 12.0, a 
bow shock will form over the forebody of the vehicle. The incoming flow first encounters a 30 (i.e., 
half apex angle) degree cone. To determine the external compression inlet performance, it is assumed 
that the surface flow conditions aft of the bow shock are uniform and equal to those entering the 
annular cowl duct. This permits the calculation of static pressure, total pressure and surface Mach 
number along the 30 degree cone (after the initial conical bow shock) as a function of flight Mach 
number. Additional flow compression caused by the parabolic primary optics surface is assumed to 
occur through a series of very weak shocks; therefore, only a small loss of total pressure occurs. The 
bow shock is always assumed to be attached to the cowl with (i.e. no spillage drag). The design 
Mach number for the inlet is 3.0, and thereafter, the shroud translates aft such that the bow shock 
is always attached to the shroud forward lip. The mass flow rate entering the vehicle is given by 
Eqn. (1). 


— PooUoo-Ac 


[III- 1) 


Poo = ambient air density 
= speed of vehicle 
Ac = capture inlet area 
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Figure 25. Scramjet configuration. 
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B. INTERNAL DIFFUSER 


In order to reduce internal cow volume, no diffuser was invoked aft of the cowl entrance 
station. As mentioned earlier, the flow conditions entering the cowl are considered to be the same as 
those after the bow shock. No normal shock is assumed to form at the cowl lip, as would be typical 
for ramjet operation. 

C. COMBUSTION CHAMBER 

Laser energy is absorbed into the flow within the combustion section, which is assumed to 
have a constant cross-sectional area. The calculations assume 85% of the incoming beam power 
appears as an increase in total temperature. The analysis computes the effect of heat addition with 
Rayleigh line equations for a constant area flow. The ratio of specific heats is set to 1.3 due to the 
elevated gas temperatures. 

Theoretical thermally-choked conditions are computed from the cowl entry Mach number. 
Laser energy then raises the total temperature of flow according to Eqn.(2). 

tie = combustion efficiency (.85) 
m == air mass flow rate through engine 
Q = laser power delivered to engine 
Cp = Specific heat at constant pressure 

Heat addition is limited such that the thermally- choked condition cannot be violated. Also the 
static temperature is not allowed to exceed 10,000 R, beyond which the ideal gas assumption is no 
longer valid. Above 10,000 R, a large fraction of the input laser energy would go into dissociation 
and ionization losses, which is undesirable. Input beam power is decreased accordingly to hold the 
temperature below these limits. 

D. PLUG NOZZLE 

In the ramjet mode, flow leaving the combustion chamber is not thermally choked; thus, a 
slight contraction in the annular duct area is needed to accelerate the flow up to Mach 1.0, prior 
to expanding over the external plug nozzle. Such nozzle area changes are generally not required for 
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scramjet operations, since the flow always leaves the combustion chamber at Mach numbers greater 
than (or equal to) one. Losses in total pressure due to this contraction are assumed negligible in 
comparison to other losses. 

The flow leaves the shroud and then expands across the bottom surface of the Lightcraft. 
The 45 degree apex half angle nozzle produces a Anal exhaust velocity determined by Eqn. (3). 

v. = {2c,n.®,[i-&) a ?*|}» (in-i) 

roc 


V t = exit exhaust velocity 
Toe = total temperature after heat addition 
r\ n a nozzle adiabatic efficiency (.95) 

Poo = ambient air pressure 
Poe a total pressure after combustion 
7 a specific heat ratio (1.3) 

The net thrust produced by the engine, neglecting the external drag on the shroud and 
pressure terms, is given by 

T*«< = mVe - mUoo {III -4) 

Tnet = net thrust 
m a mass flow rate (air) 

where the first term is referred to as gross thrust , and second is the ram drag. Vehicle drag force is 
computed using the vehicle drag coefficient, which is shown as a function of flight Mach number in 
Fig. 26. 

E. RAMJET 


A performance analysis for ramjet operation was performed for low supersonic speeds rang- 
ing from Mach 2 to Mach 6. A normal shock was assumed to be attached to the cowl lip at all 
times. Conditions across the normal shock were computed, and Rayleigh line equations to compute 
the effects of heat addition were used (as with the scramjet). 

The results indicate that the ramjet does not give positive thrust until it reaches altitudes 
and speeds where the scramjet already provides better performance. It was thus decided to abandon 
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Figure 26. Vehicle drag coefficient. 
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the ramjet mode since the seramjet performs better throughout the entire Mach 3-11 regime. The 
major reason for poor ramjet performance is intimately tied to the cowl design; no internal diffuser 
was included and large losses in total pressure were caused by the normal shock. 

F. RESULTS 

A computer code was assembled to facilitate the analysis of seramjet performance through- 
out the flight envelope. The program computed thrust, coupling coefficient, propulsive efficiency, 
laser power and overall efficiency as a function of flight Mach number, with altitude as the parametric 
variable. Graphs of these results are presented in Figs. 27 through 31. 
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Figure 27. Thrust vs. Mach number for the scramjet. 
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Power vs. Mach No. (4Gw) 
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Figure 29. Power vs. mach number for the scramjet. 
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Propulsive Efficiency 


Propulsive Efficiency vs. Mach No. (4Gw) 
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Figure 30. Propulsive efficiency vs. Mach number for the scramjet. 
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Overall Efficiency 


Overall Efficiency vs. Mach No. (4Gw) 
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Figure 31. Overall efficiency vs. Mach number for the scramjet. 
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CHAPTER IV 
MHD-FANJET MODE 


The MHD-fanjet is an electric air-turborocket which uses electromagnetic field propulsion 
to accelerate the partially ionized air heated by the bow shock formed over the forebody of the 
vehicle. The basic cycle is illustrated in Fig. 32. The laser absorption chamber procesess high 
pressure (30 to 80 atm.) liquid hydrogen propellant into a high temperature (15,000 to 20,000 K) 
ionized gas using a standing laser-supported combustion (LSC) wave. This high temperature, high 
pressure plasma is then driven via a pressure gradient through an MHD generator, which extracts 
electrical power from the plasma by applying a magnetic field perpendicular to the core flow (the 
velocity outside of the boundary layers in the generator duct) and placing electrodes along opposite 
walls of the generator duct through which the induced current can flow. The load for each pair of 
generator electrodes is a connection across the MHD accelerator duct. This current flow across the 
accelerator duct creates the applied electric field, E tn u for the MHD accelerator. A magnetic field 
is applied perpendicular to E» f p j, creating a Lorentz force which accelerates the ionized air flowing 
through the accelerator. This provides the majority of the thrust during the MHD fanjet mode, the 
additional thrust resulting from the high velocity hydrogen exhaust of the MHD generators. 

A. LASER ABSORPTION CHAMBER 

A continuous (CW) laser beam from the remote power source is focused into the absorption 
chamber to a rectangular cross-sectional area small enough (or an intensity large enough) to achieve 
electrical breakdown in the hydrogen propellant filling the chamber. The laser wavelength is chosen 
such that the propellant velocity equals the propagation velocity of the LSC wave (toward the 
laser source), producing a standing LSC wave at the inlet of the MHD generator. The absorption 
chamber/generator concept is depicted in Fig. 33. For the high aspect ratio at the inlet of the 
generator, the LSC wave will practically be a line. If this proves to be an unstable geometry, one 
possible alternative is to use a linear series of smaller circular cross-section LSC waves across the 
inlet of the generator. 

For laser powers larger than 10 kW and static pressures larger than 3 atm, the hydrogen 
LSC wave will absorb nearly all of the laser power (over 98%). However, radiation heat loss from 
the LSC wave can seriously reduce the net power going into the gas. In a numerical study by Jeng 
and Keefer [8] of an axially-symmetric hydrogen LSC wave, a wave velocity of a 100 m/s, a static 
pressure of 3 atm, and a laser power of 10 to 60 kW produced a conversion efficiency of laser power 
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Figure 32. Diagram of the MHD-fanjet cycle. 
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Figure 33. Concept for Laser-heated MHD generator. 




to net power into the flow of 65%. For the large static pressures and high velocities at the inlet 
of the MHD generator, and using highly reflective walls inside the absorption chamber, it is hoped 
that a conversion efficiency of 80% can be achieved. 

A theoretical analysis of the absorption chamber for the Apollo Lightcraft was not per- 
formed; however, A. Setagesh of Tuskegee University is investigating hydrogen LSC waves in rectan- 
gular geometries similar to the inlet of the MHD generator. The laser energy per mass of propellant 
required to raise the temperature of the hydrogen entering the absorption chamber to the inlet 
temperature of the MHD generator is currently calculated by 

q l = -L{HOi-Ho eh ) ' (rv-i) 

Vch 

where the stagnation enthalpies HO c k and HO, are calculated at the entrance conditions of the 
laser absortion chamber and MHD generator respectively, and a conversion efficiency, rich, of 100% 
is used. A more realistic value for rich will be used after an analysis of the absorption chamber has 
been completed. The liquid hydrogen, stored as para hydrogen, will be used to regeneratively cool 
the walls of the absorption chambers and MHD generators, which are assumed to support a constant 
wall temperature of 2500 K. 

B. MHD GENERATOR 


1. Geometry 

The geometry of the MHD generator is shown in Fig. 34. The cross-sectional area of the 
generator duct is rectangular, with a constant width D y and a variable height D z . The core flow 
U e is in the x-direction, the applied magnetic field B is in the z-direction, and the induced Faraday 
current density, resulting from the flow of the ionized hydrogen perpendicular to the applied 
field B, is in the negative Y-direction. The electrode walls (the two walls perpendicular to the 
current flow J v ) are segmented. The external load across each pair of electrodes is such that the 
loading parameter, K, remains constant along the generator duct. The ionized hydrogen exits into 
the ambient air aft of the vehicle contributing to the vehicle thrust. 

2. Theory 

The high pressures which exist throughout the generator allow the gas to be considered 
in local thermodynamic equilibrium (LTE); so, although the gas consists of several species, its 
physical properties may be described by a single temperature and pressure. Also, for the conditions 
encountered in the MHD generator, the gas as a whole may be treated as a continuum. Thus, the 
flow may be adequately modelled by a single set of global continuum equations consisting of the 
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conservation of mass, electrical charge, linear momentum, and energy, and a transport equation 
for the current density (usually referred to as the generalized Ohm’s law in one of its simplified 
forms). To describe the electromagnetic field interaction in the gas, the Maxwell equations for the 
curi of the magnetic intensity vector S and the curl of the electric field vector f?, along with their 
constitutive equations, are required in general. For gaseous mediums, the relative permeability p 
and permittivity e of the gas are approximately equal to their free space values of p 0 and e 0 > 

The Debye length may be interpreted as the distance by which electrons and ions will 
separate in a ionized gas due to the thermal velocity of the electrons. The Debye length (in S.I. 
units) is defined by 

ho -(?*£)* -«\A7* (TV -2) 

where T is the gas temperature and n« is the electron number density. The dimensions of the 
generators considered for the Apollo Lightcraft are on the order of centimeters (or more), while the 
maximum Debye length typically encountered is on the order of 0.1 micrometers; thus, the gas is 
essentially electrically neutral throughout the generator. 

For many MHD generators, it is often true that the average of the applied field over the cross- 
sectional area of the duct is relatively undistorted by the internally produced magnetic fields resulting 
from the electrical currents within the gas. The Maxwell equation for the magnetic induction, B, 
consistent with the assumptions of zero net charge density (and therefore the diplacement vector D 
is negligible) and a relative permeability of /*<, is given by 

$ x B = Mo/ (/V- 3) 

An order of magnitude approximation of this equation leads to a parameter describing the distortion 
of the applied magnetic field given by 

dB _ noJyD z /2 

~B B~ 

where J 9 and D z f 2 are characteristic values for the current density perpendicular to the applied 

magnetic field B and the length over which the magnetic field acts, and dB is a measure of the 

distortion or change in the magnetic field due the internal currents J y . As shown in Fig. 35, the 
parameter dB/B is approximately equal to the tangent of the angle through which the applied field 
is deflected by the internal magnetic fields. It is assumed in this analysis that for dB/B less than 
1/2, the distortion of the applied magnetic field is negligible. With this assumption, the magnetic 
field distribution becomes a known quantity and the Maxwell equation (3) is not needed. 

The transport equation for the current density may be derived by considering the first 
moment of the Boltzmann equation for each species in the gas. This leads to a highly nonlinear 


(TV -4) 
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partial differential equation for the current density vector /. However, for the MHD generator 
operating conditions, many of the terms are negligible and the transport equation for the current 
density or the generalized Ohm’s law becomes 

!=<t{$+9 xS)- ^(/x S) (IV - 5) 

where the conductivty a accounts for both electron-ion collisions and electron-neutral collisions, 
and the Hall parameter w«r« represents the ratio of the cyclotron frequency of the electrons to the 
collision frequency of the electrons. 

Also, the gas flow is assumed to be steady and will almost certainly be turbulent in any 
real application, especially for the high velocities typical of the generators considered here. The 
radiation stress tensor and internal energy, which must in general be included in the momentum 
and energy equations, can be shown to be negligible based on order of magnitude calculations for 
a black gas at typical generator temperatures (10,000 to 20,000 K in this paper). The radiative 
heat flux, however, is quite significant and will eventually be considered in the analysis; however, in 
the present analysis, radiative heat transfer between the hydrogen gas and its surroundings will be 
neglected. 

To simplify the analysis, the governing equations are integrated or averaged over the cross- 
sectional area of the generator duct. The applied magnetic field B is a function of the axial coordinate 
x only (an axially decreasing field is used) and is in the z-direction as indicated previously. The 
shear stress and heat flux are assumed uniform around the perimeter of the generator duct at a 
given x location, and the density and pressure are assumed to be functions of x only. The electrode 
walls are assumed to be infinitely segmented so that no Hall current J 3 can flow, although a Hall 
electric field E s wiil be present. Finally, making some assumptions consistent the turbulent and 
boundary layer nature of the flow, the integration of the equations over the cross-sectional area of 
the generator yields the following set of quasi-one-dimensional MHD equations: 


Continuity : pUA = m 
or 


dp ( du t dA _ q 
p 8 A 

(IV -6) 

Momentum : pUj^ + = -(»W j + JyB) 

(IV- 7) 

Energy : = ~(Qvj + JyE v ) 

(IV- 8) 


25 



(/V-9) 


Generalized Ohm's Law : J y = ^(hB - E y ) 


Non - uniformity factor : G = TT^) + (w«r«) 2 [zr(^) - l] 

(IV - 10) 

p 

Gas Law : p = < MW > 

(IV -11) 

E 

Loading parameter : K = ^ 

(IV - 12) 


Expressions for the wall shear stress r w and the convective heat flux flowing into the wall, 
q w , were obtained by approximating the flow along each wall of the generator as a steady, com- 
pressible, turbulent flat plate boundary layer flow. The friction coefficient, <7/, was calculated using 
Christoph’s and White’s inner variable method applied to a turbulent flat plate boundary layer flow 
[17], while the Stanton number was calculated using the Reynolds analogy, 

St = lc f {Pr*)W (IV - 13) 

The wall shear stress and convective heat transfer are given by 

= \pcV 2 c Cj (IV - 14) 

TT2 

q« = pU e (h e + ^.~K)St (IV - 15) 

The mole fraction for each species in the hydrogen plasma at a given temperature and 
pressure was found by applying the law of mass action to the following reactions only : 

Hi ^ 2 H 

H + H+ + e- 

With the mole fractions known, the density is calculated using a equation (11), where < MW > is 
the average molecular weight of the mixture given by 

<MW >=Y,XiMWi (IV- 16) 

i 

Although the composition of the hydrogen predicted by the above analysis neglects several important 
species and differs significantly with Patch’s data [ll], whose model considered spin-equilibrated 
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hydrogen in chemical equilibrium in the Debye - Huckel approximation and accounts for the species 
H~, H£, and values for the average molecular weight and the electron mole fraction (important 
for the calcuation of conductivity) are within 10% or better. 

The enthalpy as a function of temperature and pressure was calculated from a curve fit of 
Patch’s data, whose values are within 10% . 

The conductivity a appearing in the generalized Ohm’s law is given by the following equation 


* meC e [^n fc g fc +3.9n,( 5 ^ n .)21nAl 


(IV - 17) 


where Qk is the collision cross section of the neutral species k for momentum transfer with an elec- 
tron, C„ is the mean random velocity assuming a Maxwell-Boltzmann distribution for the electrons, 
and A is the ratio of the Debye shielding length to the impact parameter for 90° Coulomb scattering. 

The viscocity of the gas was approximated by a curve fit of data obtained by Yost [19]. 

3. Generator Design 

Equations (6) - (8) represent a set of coupled, nonlinear, first order differential equations, 
which can be expressed, at least implicitly, in terms of 7, P, 8, D z , D y and x. Typically, one 
would consider x as the independent variable for some prescribed area distribution, A = D v D z (x). 
However, in order to design a generator for the Apollo Lightcraft, it is advantageous to leave the 
area distribution variable (but with a specified shape, ie., a rectangular cross section) and specify 
some other condition to close the set of equations, since the volume of the generator can greatly 
effect its efficiency and power output, while its surface area will effect the heat and friction losses at 
the walls of the generator. The condition imposed in this analysis is that of a constant velocity, 8, 
throughout the generator. This simplifies the solution of the set of differential equations describing 
the generator, but it also removes the threat of flow separation, a problem which may occur in a 
flow with a decreasing velocity. 

The solution to this set of differential equations was performed by approximating each 
differential in Eqns. (6) • (8) by a forward difference: 


<*/(*) /(* + As) - /(ar) 


dx 


Ax 


(j TV - 18) 


The differential of the enthalpy can be written as 

ih^^'ST + ^iP^CrdT + CTiP (/V - 19) 


Cp and Ct were also approximated by forward differences. Using T as the marching variable, Eqns. 
(6) - (8) can be solved to yield: 


Ax = 


-CpAT 

[JyEy +?«£)£ - Cr{JyB + r«,£) 


(IV - 20) 
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{IV- 21) 


APm-{J,B + T W £)Ax 

For a given set of initial conditions, we calculate all the necessary variables, let T increase 
by some AT (ie., we march in T), and calculate the change in x and P using Eqns. (20) and (21). 
p and D z are updated according to 

p(x + Ax) = p{T{x) + AT, P(r) + AP) {IV - 22) 

D> -ss- 

In the above equations and in any futher references, the bar superscript over averaged quantities 
will be left out for convenience. Also, at each step, the electrical work extracted per unit mass, Aw, 
and the convective heat loss to the walls per unit mass, A(QT), are calculated by 

4(<3r)= s „(i)£jjj£ (iv -u) 

At o = h(x) - h{x + Ax) - A(QT) {TV - 25) 

The initial conditions for the generator consist of the entrance temperature, pressure, Mach number, 
and magnetic field, the loading parameter, and the mass flow rate. The parameter dB/B and the 
velocity boundary layer thickness (assumed the same for each wall) are also calculated at each step 
through the generator. 

Preliminary trajectory calculations indicated that a specific electric work of about 700 M J /kg 
would be necessary if the Apollo Lightcraft was to achieve low Earth orbit and deorbit using the 
design point of 330 kg of hydrogen propellant or less. Also, the size of the generator was limited to a 
maximaum axial length of L=0.5 m to meet with the design point volume for the Apollo Lightcraft. 
As depicted in Fig. 39 (and in the Appendix), the vehicle will contain twelve (12) MHD generators 
which extend out during the MHD fanjet mode. Since the exterior wall in which the generators 
are contained must be parallel with the shroud, the ratio L /D« (where D« is the exit height of the 
generator) must be greater than 1.5. In order to maintain the generator walls at a temperature of 
2500 K, using only the 330 kg of liquid para hydrogen to cool the walls of the generators and the 
laser absorption chambers, the total convective heat/mass into the generator walls, QT, must be 
less than some maximum heat/mass, It is also desirable to have a high conversion efficiency 

of laser energy into total generated electrical energy, w(tot)/Q^. 

The design of the generator for the Apollo Lightcraft was accomplished by structuring the 
above into the following optimization problem : 
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MAXIMIZE the objective function f = w(tot) Subject to : 


w > 70QMJfkg 
L < 0.5m 
L/D e > 2 
QT < QT m tx 
w(tot)/Qi > 0.5 

Pe > OMQSbart 


plus 

dB/B < 0.5 
2 6/D m < 1 
2S/D y < 1 

which must be satisfied throughout the duct. 

It should be noted that large ratios of LfD e help reduce losses at the generator exit due to 
eddy current loops. The constraint on the exit pressure is necessary for a generator with subsonic 
flow at the exit to function at a minimum altitude of 100,000ft. If the flow becomes supersonic within 
the generator and the exit pressure is less than the ambient pressure, a shock wave may form in the 
generator, which would presumably degrade generator performance. The last three constraints are 
imposed to ensure the validity of the theoretical model described previously for the analysis of the 
generator. 

A generalized reduced gradient (GRG) code (OPT 3.1, Gabriel and Beltracchi (5]) was used 
to solve the above optimization problem. The design was attemped for two different expressions 
for QT m tx. In what will be referred to as the HIGH case, QTmtx was set to a constant value of 
28 MJ/kg. We can place a maximum value on QTmtx for a given storage temperature of 20 K for 
the liquid hydrogen by computing the heat/mass that could be absorbed by the hydrogen if it was 
heated to 2500 K; assuming a constant C v of 15,000 J/kg/K, QTmtx would be approximately 37 
MJ/kg. The above value of 28 MJ/kg was used since the laser absorption chambers must also be 
cooled. 

In what will be referred to as the LOW case, the value for QTmtx was computed assuming the 
cold hydrogen flowed through a simple cooling duct 1 mm in height along each wall of the generator 
and assuming a certain amount of heat loss in the laser absorption chamber estimated from a mean 
beam length calculation. In this case, QTmtx depended on the geometry of the generator and the 
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temperature and pressure of the hydrogen at the generator inlet, and thus a separate subroutine 
for calculating QT mxx was incorporated into the previously described generator code. Typically, 
QT mM was on the order of 2 to 6 MJ/kg. This case is considered a lower bound on QT max since 
transpiration cooling or other advanced heat transfer methods would allow higher values for QT max . 

4. Results 

The resulting designs for the HIGH and LOW generators are described in Fig.'s 36 and 37. 
A feasible design point satisfying w(tot)=700 MJ/kg could not be found for the other constraints 
imposed on the Apollo Lightcraft design; thus, the optimization problem was simply to minimize 
the constraint violations. In both the HIGH and LOW cases, the constraints on the generator 
length L and total heat loss to the walls QT are active (ie., essentially equal to zero), while all 
other constraints (except w(tot) > 700 MJ/kg) were satisfied. The HIGH generator nearly satisfies 
the specific work goal of 700 MJ/kg, producing 652.5 MJ/kg; as will be discussed in the trajectory 
analysis (Chapter VI), this generator provides sufficient work/mass for the Lightcraft to achieve low 
Earth orbit and deorbit using no more hydrogen propellant than the design point of 330 kg. The 
LOW generator can only produce a specific work of 362 MJ/kg in order to satisfy the tighter heat 
loss constraint; this would require over 330 kg of hydrogen to achieve low Earth orbit and deorbit 
(discussed further in Chapter VI). The two designs attempt to provide an upper and lower bound on 
the generator performance posable for the Apollo Lightcraft: 652.5 MJ/kg and 362 MJ/kg electric 
work/mass respectively. 

5. Generator Thrust 

Consider the single generator shown in Fig. 38. Ignoring the contribution to thrust from 
the expansion of the exhausted hydrogen over the plug nozzle (to ambient pressure), the generator 
thrust was calculated as 

Tgen = [mu + (P t - Pa) A t ] cos /? {IV- 26) 


/ d = 90°-(45 <, + tf) 


9 = tan -1 ( 


Dt-Di 

21 


) 


G. MHD ACCELERATOR 

The accelerator duct is an annular channel formed between the shroud and the exterior walls 
of the MHD generators, which are located symmetrically about the base of the vehicle, as shown in 
Figs. 39 and 40. In the MHD-fanjet mode, the shroud will be in the most rear-ward position so 
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Inlet conditions and generator parameters : 


T = 20 r 000 K 

m = 

0.4850 

kg/s 


p = 80.000 bars 

K = 

0.9898 



B = 

4.740 tesla 

Dy = 

33.0000 

cm 


M = 

0.534 

u = 

8104.40 

m/s 


X(cm) 

Dz/2 (cm) 

T/1 000 

PG ( bar s ) 

QT(MJ/kg ) W(MJ/kg) 

0.00 

0.24 

20.00 

80.00 

0.00 

0.00 

4.43 

0.31 

19.14 

58.28 

7.90 

63.98 

8.76 

0.41 

18.28 

41.38 

13.48 

129.23 

12.90 

0.56 

17.42 

28.60 

17.34 

195.38 

16.88 

0.78 

16.55 

19.22 

20.01 

262.11 

20.82 

1.12 

15.69 

12.54 

21 .92 

329.05 

24.92 

1 . 64 

14.83 

7.92 

23.37 

395.88 

29.44 

2.48 

13.97 

4.84 

24.55 

462.22 

34.75 

3.86 

13.11 

2.86 

25.60 

527.58 

41.21 

6.18 

12.25 

1 .63 

26.59 

591 .31 

49.10 

10.16 

11.38 

0.90 

27.59 

652.51 

Laser 

Energy required, Q, 


= 1 . 

2364 GJ/kg 

Maximum value of 

dB/B 


= 0 . 

0364 

Conversion efficiency, w(tot)/Q u % 

= 52 

.78 % 


Figure 36 . The HIGH case MHD generator. 




Inlet conditions and generator parameters : 


T = 

14,916 K 

m = 

0.6974 kg/s 


P = 

30.346 bars 

K = 

0.9900 



B = 

9.930 tesla 

Dy = 

50.6900 cm 



M = 

0.353 

u = 

4311.30 m/s 



X(cm) Dz/2(cm) 

T/ 1 000 

P( bars ) QT(MJ/kg) 

W(MJ/kg) 

0.00 

0.71 

14.92 

30.35 

0.00 

0.00 

2.07 

0.95 

14.22 

21.43 

0.57 

44.79 

3.95 

1 .28 

13.53 

14.93 

0.98 

88.63 

5.74 

1 .74 

12.84 

10.28 

1 .29 

131.22 

7.52 

2.39 

12.15 

7.01 

1 .55 

172.21 

9.42 

3.30 

11 .46 

4.75 

1 .78 

211.19 

11.68 

4.55 

10.76 

3.21 

2.02 

247.74 

14.68 

6.23 

10.07 

2.18 

2.30 

281.45 

19.36 

8.43 

9.38 

1 .49 

2.69 

31 1.96 

28.28 

n . 1 8 

8.69 

1 .04 

3.39 

339.01 

50.26 

14.42 

8.00 

0.74 

4.97 

362. 19 

Laser 

Energy required , Q, 

= 

0.7219 

GJ/kg 

Maximum value of 

dB/B 

zz 

0.0045 


Conversion efficiency, w(tot)/Q u % = 

50.17 % 



Figure 37. The LOW case MHD generator. 




that the forebody bow shock is attached at all times. For Mach numbers above 12, the position of 
the shock wave does not change significantly. The principle of operation for the MHD accelerator is 
unique in that energy is added to the flow through an electromagnetic field and converted directly 
to kinetic energy. For vehicles travelling in the hypersonic regime, the formation of a strong bow 
shock can result in static temperatures downstream of the shock on the order of several thousand 
degrees Kelvin. As a result, the air passing through this region is partially ionized and effective 
electrical conductivity is increased. This conductivity is then in general a function of altitude and 
Mach number, as well as bluntness of the vehicle forebody. 

The flight envelope currently envisioned for the MHD fanjet mode begins roughly at Mach 
12 at 100,000 ft. This approaches the limit of operation of the Scamjet. Obtaining analytical 
performance data for this engine involved an investigation of the flow over the conical forebody, 
and the subsequent acceleration in the MHD channel itself. Operation at the relevent hypersonic 
velocities required the use of real gas data both through the inlet and in the MHD accelerator. 

1. Theory 

The assumptions made for the hydrogen flow through the MHD generator also apply to 
the air flow through the MHD accelerator; thus, the quasi-one-dimensionai equations (6) - (12) can 
be used. Some father simplifications are made for the MHD accelerator analysis. The shear stress 
and heat transfer (both conductive and radiative) at the wall are neglected, as well the Hail term, 
JxS, in Ohm’s law (5). Also, the differential of the enthalpy is represented by dh = C r dT, where 
C p is considered a function of temperature and pressure. For the accelerator, the applied magnetic 
field and the cross-sectional area are constant throughout the accelerator channel. With these 
conditions and assumptions, the quasi-one-dimensionai MHD equations for the MHD accelerator 
analysis become: 


Continuity : ^ = 0 

P 

(IV- 27) 

Momentum : pu^- + — = J y B 
dx dx 9 

(IV - 28 ) 

Energy : puC F ^ + pu 2 ^ = J y E y 

(IV - 29) 

Ohm's Law : J y — a(E y - UB) 

(IV - 30) 

Equation of State : R = < ^ > ; p = ^ 

(TV - 31) 


The parameter K is still defined as E y / uB, but for the case of the accelerator is greater than 1. 

2. Accelerator Model 

Re-arranging Eqns. (27) - (31) in terms of differentials for velocity, temperature, pressure, 
and density, using Ohm’s law and the loading parameter to eliminate J y and E y , and defining the 
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parameter Z as the mass flow rate per unit area (constant), gives the following set of relations: 


«•*{*- 0(1 -») 

(IV- 32) 

r _ <tv?B 2 K(K - l)Ax uAu 
ZCp Cp 

(IV - 33) 

AP = ^jp[AT - — Auj 

(IV - 34) 

A 1 lAO P&T, 

Ap= wl AP ~—\ 

(IV -35) 


These equations are readily integrated for prescribed values of B and K, initial values of u, 
T and P, and a step size Ax. An inlet analysis provided the required initial velocity and thermody- 
namic state variables, and a tabulation of the electrical conductivity, specific heat, and gas constant 
provided the necessary thermophysical properties. 

Data for C p , a, and R for air at elevated temperatures over a wide range of pressures was 
obtained from data calculated by Yos (18|. While not tabulated explicitly, it was possible to calculate 
the value of the gas constant R from the given data for temperature, pressure, and density. These 
tables were then made available to a subroutine which would determine the required properties for 
the given temperature and pressure at each step of the solution. Fig. 41 shows a typical profile 
solution in non-dimensional form. 

3. Isothermal Solution 

Heretofore, nothing has been mentioned regarding the selection of the loading parameter K, 
except that it is prescribed for a given solution. The applied magnetic field is held constant along the 
length of the accelerator for a particular solution. If the loading parameter is also held constant, the 
result is an electric field profile increasing in proportion to the velocity. More generally, it is possible 
to alter the shape of the velocity profile by varying the applied electric field. It should therefore be 
possible to select the electric field profile to optimize the velocity increase through the accelerator 
for a given magnetic field. This analysis was done by Drake [4], who uses techniques of variational 
calculus to show that the solution corresponding to minimum channel length is equivalent to the 
solution for maximum thrust per unit exit area. More importantly, the optimal solution is found to 
be very close to that for a duct of constant cross-sectional area, which is the case of interest in the 
present analysis. 

Ideally, one would want all the applied electrical energy to be converted into kinetic energy. 
The implication is that the flow will be isothermal since no energy goes into heating the gas. This 
requirement then provides the means of selecting K (and ultimately E v ) as a function of x by 
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requiring the flow to be isothermal. This relation is readily obtained by setting AT equal to zero in 
Eqn. (33), eliminating Au with Eqn. (32) and solving for K. This gives: 

K TTB[ (/V - 36) 

which is the expression for K, and hence to insure the temperature remains constant throughout 
the accelerator. The electric field profile representative of the flow solution in Fig. 41 can be seen 
in Fig. 42. 

As a practical matter, obtaining such a distribution along a constant area accelerator as- 
sumes the use of segmented electrodes, each at a different potential. Such segmenting, while perhaps 
difficult from the standpoint of insulating adjacent segments, is nevertheless a desirable feature to 
reduce hall currents. 

4. Inlet Analysis 

As was mentioned previously, before the numerical solution can be initiated, it is necessary to 
prescribe a set of initial conditions. These are obtained from analysis of the flow through the conical 
bow shock. For the range of Mach numbers under consideration (> 12), the static temperatuures 
even behind an oblique shock preclude the use of ideal gas relations for determining the temperature 
and pressure downstream of the shock. Oblique shock data for real gases can be found in the 
literature [16] and curve fits based on this data were incorporated into the engine model to calculate 
temperature and pressure rises across the shock. With static temperature and pressure known it 
is possible to find the density using the equation of state and real gas properties for the given 
conditions. Finally the velocity is determined from the equation of continuity at the accelerator 
inlet. 

The flow field about the conical inlet will be self similar with properties constant along 
conical surfaces swept out by rays emanating from the vertex. For the 1-D analysis, properties were 
assumed constant across the inlet plane of the accelerator. This is a reasonable assumption provided 
the channel height is small compared to the cone length. 

D. MHD-FANJET PERFORMANCE 

With the inlet and accelerator models established, it is possible to calculate the thrust, 
specific impulse, and coupling coefficient for the MHD-fanjet. 

f = m*(Ue ~ Vi) + (P e - Pi)A F + (12 r gtn )THROT 
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(IV - 39) 


T 

Pl 

mj — (12 m)THROT 

THROT- Tf Mca p acity 

where mj is the total hydrogen mass flow rate for the MHD generators, r jen is the thrust contribution 
from a single generator exhaust, and THROT is the throttle setting for the MHD generators. It 
is assumed for a given generator that the specific work, w(tot), and the conversion effieciency, 
w(tot )IQl, remain constant as the mass flow rate through the generator is decreased somewhat 
about its design value The electrical power required by the accelerator is found by integrating the 
electrical power input per unit volume over the volume of the accelerator 

P F = J J J J y E y d(V ol) = B*A f J' <jv?K(K - l)dx (IV - 40) 

Curves for thrust, specific impulse, and coupling coefficient are shown in Fig. 43 over a 
range of altitudes and Mach numbers, for the case of B=2 Tesla and an available laser power of 
10GW (the maximum power level for the Lightcraft). The curves are representitive of the typical 
behavior of the MHD-fanjet. A characteristic of MHD-fanjet mode is that fuel consumption (in this 
case, the hydrogen in the MHD generators) is constant for a given power level regardless of altitude 
or Mach number. As a result, the speciflc impulse and coupling coefficient are proportional to the 
thrust over the entire flight envelope, and it is possible to represent these three parameters with a 
single family of curves. The speciflc impulse of the MHD-fanjet is generally four to five times greater 
than that of the laser rocket (I sp = 2000 s). 

It is important to note that for the isothermal operation of the accelerator, the electrical 
conductivity (which is a strong function of temperature and a weaker function of pressure) is very 
nearly at the value corresponding to the static temperature at the inlet of the accelerator. Through- 
out most of the flight regime, the electrical conductivity coiTresponding to the inlet temperature is 
not high enough for the MHD accelerator produce positive thrust. For this reason, it is necessary 
to enhance the conductivity through some additional means, such as microwave or radio frequency 
exitation, electron beams, or simple chemical seeding. The advantages and disadvantages of each of 
these alternatives is currently being investigated. In this model, the electrical conductivity, (Jtot, 
required to absorb the power available from the MHD generators was computed by 

&TOT — (T + (T g 

where a t is a constant representing the amount by which the actual conductivity a must be increased 
such that P F (found by Eqn. (40), with otot replacing <r) is equal to the available power produced 
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by the MHD generator*, m/w(tot); <r, i* found iteratively. The maximum value of <r # , at a given 
Mach number and altitude, occur* at the maximum laser power of 10GW. <r t verses flight Mach 
number at various magnetic field values is plotted in Fig. 44 for a laser power of 10GW at an 
altitude of 170,000ft. 
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CHAPTER V 
ROCKET MODE 


At altitudes greater than 200,000 ft, the MHD-fanjet begins to function essentially as a 
rocket since the thrust from the MHD generator hydrogen exhaust begins constitute a large fraction 
of the total MHD-fanjet thrust, due to the decreasing ambient density. Consequently, at an altitude 
of 260,000ft, the applied magnetic fields for the MHD generators and accelerator are shut off, and the 
Lightcraft transitions into a pure rocket mode - the twelve MHD generators become laser-heated 
hydrogen rockets. The laser absorption chamber and MHD generator duct form a converging- 
diverging nozzle in which the throat area is the inlet of the MHD generator. In the trajectories 
currently used, the accelerators and generators are shut off at 260,000ft. From here, the Lightcraft 
coasts out to about 600,000ft, at which point, a rocket burn is used for orbital insertion. However, 
in future studies, the use of an additional rocket burn at the end of the MHD-fanjet mode may be 
investigated as a possible strategy for minimizing total energy consumption. 

A. ROCKET PERFORMANCE 

For a given MHD generator (ie, the HIGH or LOW case generator), the geometry of the 
generator duct and laser absorption chamber are fixed. However, since the laser absorption cham- 
ber temperature and pressure can be varied about their MHD-fanjet values (by altering the laser 
energy into the chamber and the pumping pressure of the hydrogen propellant), and the absorption 
chamber/generator duct constitutes a converging-diverging nozzle, the MHD generators can be used 
as rocket nozzles, assuming the back pressure (ambient pressure) is low enough to allow supersonic 
flow at the exit of the generator duct. As noted, the minimum altitude at which the rocket mode 
will be used is 260,000ft, which corresponds to a maximum back pressure of 1.74 x 10“ 5 atm; thus, 
for entrance stagnation pressures on the order of 1 atm (or more), the flow will be supersonic and 
underexpanded at the generator exit, with further expansion occuring over the plug nozzle. 

A detailed analysis of the hydrogen flow through the generator duct in the rocket mode 
was not performed, although such an analysis will soon be carried out using the model developed 
for the MHD generator in the MHD-fanjet mode, with the electromagnetic terms removed and the 
constant velocity condition replaced by a prescribed area distribution for the generator duct. In 
order to carry out the trajectory analysis for the rocket mode, some engineering assumptions were 
made. The thermal efficiency of the generator duct operating as a rocket was assumed to be 66 % 
(reasonable for the high operating temperatures) and the thrust from each rocket was defined as the 
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product of the hydrogen mass flow rate times the exit velocity (the expansion over the plug nozzle 
was ignored). Defining I» v as the specific impulse and C as the coupling coefficient, the thermal 
efficiency, can be written as t jt* = where g is the gravitational acceleration. We can 

generate very high specific impulses for the rocket since high temperatures (15,000 to 20,000 K) 
are easily obtained using standing LSC waves in the absorption chamber; however, in order to also 
obtain a reasonable coupling coefficient of 67 N/MW, a specific impulse of 2000 sec was decided on 
as the design point for the laser-heated rockets. 

The desired specific impulse can be generated for a given generator duct by providing the 
proper stagnation temperature and pressure at the absorption chamber inlet. Each rocket must 
also produce a minimum thrust of approximately 18.6 kN to satisfy certain requirements for the 
optimization of the trajectory (assuming all 12 rockets are operating at the same capacity). The 
mass flow rate for each rocket must therefore be greater than or equal to 0.948 kg/s for a specific 
impulse of 2000 sec. For a given TO (Stagnation temperature) and PO (Stagnation pressure), the 
mass flow rate per unit area necessary for a choked flow is determined, and for a given generator 
inlet area, Ai = DiD v (ie, the throat area of the rocket nozzle), the mass flow rate is determined. 
Using rocket performance data for high temperature hydrogen generated by Patch [?], it was possible 
to check the feasibility of achieving a specific impulse of 2000 sec and a minumum thrust of 18.6 
kN for the HIGH and LOW case generators. For the HIGH generator, the exit to throat area ratio 
is 42. For T0=8000 K and PO=60 atm, linear interpolation of Patch's data yields: /« p = 2040s, 
an exit pressure Pe of .33 atm, and a mass flow rate m of 0.898 kg/s. Using the same conditions 
for the LOW generator (which has a exit to throat area ratio of 20.4) yields : = 1955s, Pe = 

0.52 atm, and m = 4.09 kg/s. Of course, these calculations do not account for the effect of the wall 
shape between the throat and exit of the nozzle, nor the heat and friction losses at the walls. They 
do indicate, however, that the numbers used for the specific impulse and coupling coefficient in the 
trajectory analysis for the rocket mode are reasonable. 
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CHAPTER VI 
TRAJECTORY ANALYSIS 


The trajectory analysis of any launch vehicle becomes a critical step in the overall system 
integration process. Many important engine/vehicle related characteristics of the craft must come 
together for the find product - performance. Launch vehicle performance is typically measured in 
terms of payload capability, which for the Apollo Lightcraft is the five person crew (plus space suits) 
which totals 500 kilograms. This vehicle must be capable of attaining low Earth orbit with the 
available amount of propellant while minimizing the total laser energy consumed along the insertion 
trajectory. 

The trajectory was evaluated using a computer tool called SORT (Simulation and Optimiza- 
tion of Rocket Trajectories) which was written by McDonnell Douglas Astronautics Co. for NASA 
to design space shuttle trajectories (2j. The code is sufficiently general such that it can analyze any 
trajectory (even those for other planets), and model all important environmental parameters that 
affect the vehicle dynamics. Engine performance, vehicle aerodynamics, guidance algorithms, and 
mass histories interact with with atmosphere and gravity models. 

The capability of the SORT program is quite extensive in its use of these sophisticated 
vehicle and environmental models. The program can iterate on trajectory parameters to optimize 
performance, achieve a desired criteria, or constrain the solution to avoid some specified limit. These 
capabilities make SORT one of the most powerful trajectory analysis tools available in the world. 

Even with all of the generality built into SORT certain modifications were required for the 
Apollo Lightcraft vehicle. The most significant software modification involved the unique energy 
source, a laser. A new vehicle steering option was encoded so that the Lightcraft could always point 
at the laser power relay satellite which initially passes overhead (See Fig. 45). Since the use of laser 
propulsion was not anticipated by the authors of SORT, software changes were also needed in the 
engine model in order to evaluate the total amount of laser energy expended along the trajectory. 
After these modifications to the code were included, SORT was capable of modeling Lightcraft 
performance to a high degree of accuracy. 

A. VEHICLE AND ENVIRONMENTAL MODELS 

1. Aerodynamics 

For this preliminary performance analysis vehicle drag was the only aerodynamic force (or 
moment) that was modeled. Since the Lightcraft is axisymmetric, it produces no lift at zero angle- 
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of-attack. During flight the vehicle angle-of-attack stays fairly small due to the restrictive steering 
requirements and high thrust to weight ratios. It is suspected that aerodynamic lift can possibly be 
used to benefit Lightcr&ft performance along the trajectory. 

The vehicle drag profile is shown in Fig. 46 The very low drag coefficients indicated for 
most Mach numbers is a result of the vehicle configuration, which effectively is a flying engine. In 
the analysis of hypersonic air breathing launch vehicles, it is not always immediately evident as to 
what should be classified as “engine " versus “airframe ” - especially when large portions of the 
vehicle forebody provide the hypersonic inlet function. Hence careful inventory must be kept to 
avoid penalizing the engine for what might be more appropriately accounted against the the vehicle 
airframe. 

For example, the large form drag normally associated with “blunt " supersonic bodies is 
treated as a loss term in in the engine model. Almost all of the Lightcraft forebody is the supersonic 
inlet and the pressure on this surface (as well as the external nozzle) are modeled in the engine 
calculations. The shroud is the only remaining aerodynamic interaction with the freestream. The 
drag produced by this surface will vary considerably depending on the deflection or curvature of 
the shroud. The drag coefficient values past Mach 3, shown in Fig. 46, correspond to 15 degree 
deflections of the leading and trailing edges (see Fig. 47). The large drag increase from Mach .9 to 
Mach 3.0 is an artifact of “spillage ” drag. The inlet does not capture all of the of the shocked air 
until the vehicle reaches Mach 3, the design Mach number of the Apollo Lightcraft inlet. Hence, 
below Mach 3, the air which is slowed by the shock, but not ingested by the engine is attributed to 
vehicular drag. The shroud is designed to translate fore and aft, such that all of the shocked air is 
swallowed at flight speeds beyond Mach 3. 

2. Engine 

Detailed discussions for each engine model are presented in earlier sections of this report. 
The flight performance for each engine was calculated for a wide range of altitudes, Mach numbers, 
and power settings. These values were then loaded into trivariant tables accessible to the program. 
The SORT program interpolates this data to find the thrust, specific impulse and power consumption 
at any point along the trajectory. The engine model also calculates the cumulative laser energy used 
throughout the launch trajectory, which also includes atmospheric attenuation of the laser beam. 

3. Steering 

The Apollo Lightcraft must point at the laser relay satellite at all times during powered 
flight. Two laser power source locations were assumed for most of the analyses which have been 
carried out to date. One location is directly overhead in a very high orbit, where no significant 
movement occurs with respect to the earth’s surface during the boost; the other location considers 
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Figure 47. Lightcraft cowl profile. 


MASS BREAKDOWN 

FOR "APOLLO 

LIGHTCRAFT" 

COMPONENT 

MASS (KG) 

WEIGHT (LBS) 

PAYLOAD (5 PERSON CREW) 

500. 

1100. 

STRUCTURE (+ PRIMARY OPTICS) 

2550. 

5620. 

THERMAL PROTECTION SYSTEM 

450. 

990. 

SHROUD (+ SECONDARY OPTICS) 

180. 

400. 

SHROUD NOZZLE (VARIABLE AREA) 

180. 

400. 

LASER-HEATED ROCKETS 

400. 

880. 

MHD GENERATORS 

200. 

440. 

MHD AIR ACCELERATOR 

200. 

440. 

LANDING GEAR (RETRACTABLE) 

170. 

380. 

6-AXIS CONTROL ROCKETS 

100. 

220. 

CHUTE AND FLOTATION 

90. 

200. 

AVIONICS, ENVIRON. CONTROL, 
AND PERSONAL PROVISION 

120. 

260. 

LH 2 TANK (INSULATED) 

80. 

180. 

LH 2 PROPELLANT 

330. 

740. 

TOTAL 

5550. KG 

12,250. LBS 


Figure 48. Mass Breakdown for the Apollo Lightcraft. 
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a second relay satellite passing overhead in a low orbit (100-200 NM). Nearly all of the acceleration 
takes place as the vehicle chases the low altitude relay satellite. Currently, the steering model works 
only for satellites passing directly over the launch site, and further improvements should permit the 
inclusion of relay satellites with a more general ground track. 

4. Mass Properties 

The preliminary weight breakdown for the Apollo Lightcraft can be found in Fig. 48. The 
program can model variations in the center-of-gravity (CG) location, but the present analysis placed 
the vehicle CG on the centerline (along with the engine) to simplify the thrust-pointing requirements. 

5. Gravity 

A complete Earth gravity model exists in the SORT code which includes J2, J3, and J4 
effects. The program continuously calculates the orbital parameters throughout the boost trajectory. 

6. Atmosphere 

The atmospheric data loaded into SORT was the ARDC Model Atmosphere [6]. Although 
wind can be a significant aerodynamic perturbation none was included at this stage of the analysis. 

B. TRAJECTORY SIMULATION 

The SORT program treats the trajectory as a specified sequence of maneuvers or phases. 
Each event is defined by the user at appropriate places in the simulation so that new inputs can be 
assembled into the trajectory. The launch vehicle flies from event to event until, finally, it obtains 
the orbital parameters necessary to achieve the desired orbit altitude. The craft then coasts until 
it reaches apogee, and performs a final orbit-insertion bum. The SORT program has iteration 
capability that can be used to optimize total propellant weight (or laser energy) used during ascent. 
A typical trajectory sequence is described below. 

1. Liftoff in ERH thruster mode 

The vehicle climbs vertically (using the high orbit satellite) to a specified altitude. This 
“pop-up " maneuver lifts the vehicle out of the dense atmosphere before attempting the acceleration 
run. Thereby, the high drag losses and structural loads associated with very large dynamic pressures 
are avoided. Simulations prove that the pop-up maneuver is more advantageous than accelerating 
“off the deck ”, Since the ERH thruster mode is completely air-breathing, the manuever does not 
consume any LH? propellant and only a small penalty in laser energy cost is incurred - which is 
paid back later in the trajectory, because the aerodynamic forces (and hence, required laser energy) 
are reduced. 
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2. Pitchover (ERH thruster) 

At the specified (or optimized) altitude the vehicle pitches over to begin receiving power 
from the low altitude relay satellite. The exact zenith angle at which to engage the relay satellite 
is determined iteratively to achieve maximum performance. If the satellite is engaged when it is 
too high overhead, the vehicle will leave the atmosphere before high velocities are reached, which 
necessitates more use of the rocket engine. If engaged too late (i.e., when the satellite is too far 
from vertical the satellite will disappear over the horizon before orbital speeds are attained; or the 
trajectory will not climb high enough to prevent increased aerodynamic drag which is associated 
with additional expenditures of fuel or laser energy. 

This maneuver effectively initiates the acceleration portion of the flight. As the vehicle 
pitches over, the pulse repetition frequency of the ERH engine is increased to accelerate the craft 
through the transonic region towards Mach 3. As the ram drag increases net effective thrust of the 
ERH thruster decays. 

3. Scramjet Mode 

To increase vehicle thrust the next air-breathing engine mode begins and the ERH thruster 
is turned off. Currently the scramjet mode is used to accelerate the vehicle from about Mach 3 
to 11. The scramjet portion is flown using a fixed (or optimized) maximum laser power level. At 
moderate altitudes, the scramjet can handle relatively high power inputs because of large mass flow 
rates and “low " inlet temperatures. However all of the 10 GW (max laser power limit) is not 
needed. Three Gigawatts of laser power can accelerate the vehicle at up to 6 g’s (even with vehicle 
drag forces included). If the craft accelerates too quickly, the aerodynamic forces build up, and as 
a result acceleration then declines while structural loads increase. At high Mach numbers, or as the 
atmosphere thins out at higher altitudes the engine can only add small amounts of energy to the 
air flow without exceeding the temperature limits of the engine (set at 10000°R). As thrust decays 
near Mach 11 and vehicle acceleration approaches “break even " (1 G), the MHD fanjet will start. 

4. MHD Fanjet Mode 

When the scramjet fades out and acceleration drops below 40 ft/ sec 1 the fanjet is started. 
By waiting as long as possible to switch to fanjet it is possible to conserve onboard LH<i propellant. 
The fanjet mode is used to accelerate the vehicle right up to orbital insertion speeds. As long as 
it captures significant mass flow rates of air to accelerate, the MHD fanjet has specific impulses far 
above the rocket mode. Hence, the rocket is used only for final “burns " in orbit. 

The final velocity desired at fanjet termination is a function of altitude, flight path angle 
and the drag experienced as the vehicle coasts up to the desired altitude (e.g. 100 NM) to make the 
final orbit circularization burn. There is no easy analytical solution to this targeting problem so the 
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program iterates to find the fanjet cutoff velocity. 

5. Coast to Apogee 

This is the engineless phase where the vehicle coasts up to apogee. After the fanjet cuts 
off the drag coefficient is increased to .3 since the forebody drag increases significantly when the 
shocked air is not re-accelerated within the cowl. 

6. Final Insertion 

The rocket is started as the vehicle reaches apogee. This bum is terminated as the inertial 
velocity reaches the necessary value for circular orbit at that altitude. 

Many trajectories are run in any optimizing attempt. Typically the performance indicator 
propellant weight was optimized, although total laser energy was run also. Analysis in the future will 
focus on the minimization of laser energy. The key parameters which were varied in the attempt to 
optimize weight were the initial position of the relay satellite, power available to the engines (throttle 
setting), and altitude for the pitchover maneuver. Trajectory optimization is a never- ending task 
(e.g. engineers are still tweaking Space Shuttle trajectories to discover more payload capability), 
but the results so far give a baseline indication of the Apollo Lightcraft performance potential. 

C. RESULTS 

The trajectory analyses to date using SORT have just barely scratched the surface of the 
comprehensive performance optimization problem. However these simulations do give highly accu- 
rate answers for the selected models and input data. Assuming the vehicle models (engine, drag, 
etc.) are fairly accurate, the resultant performance reported here would certainly be attainable. 
Although not optimal, the current results indicate the types of problems to be encountered, and 
present conservative performance capabilities (within the accuracies of the models and assumptions). 
Further, more detailed analyses cam be expected to improve upon these results. 

The earliest trajectory concepts for the Apollo Lightcraft focussed on maximum thrust 
takeoffs, reasoning that high accelerations would give the best performance. That approach does 
not bear fruit for this advanced launch vehicle which has incredible thrust-to-weight capability. The 
“max-G " type of launch quickly builds up very high aerodynamic loads ( dynamic pressures 1 10000 
psf), and the drag rises to meet the thrust. Engine failure at these dynamic pressures could result 
in decelerations as large as 50 G’s or more! Alternatively, low acceleration takeoffs had problems 
with the satellite disappearing over the horizon before the insertion was complete. 

As the SORT program became available a much improved analysis of the trajectory began. 
It became obvious that the problem of the high drag transonic region could be eased by flying with 
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lower dynamic pressures. This eventually led to the origination of the pop-up maneuver concept, 
wherein the vehicle climbs out subsonically to a high altitude before attempting to accelerate through 
this drag barrier. 

When using a single relay satellite for a slow climb-out the vehicle would gradually pitch 
over to follow the satellite. As this occurrs, more thrust would be needed to maintain at least one 
(1) G in the vertical direction. To solve this problem, it was suggested that a very high altitude 
satellite (e.g. geosynchronous) would move slowly (or not at all), allowing the vehicle to climb almost 
vertically. This option also required very low laser power requirements, which opened the possibility 
of using microwave power for the vertical climb. Since microwaves are not significantly attenuated 
by clouds, this would allow all-weather operation. Even on a cloudy day one could fly to orbit in a 
lightcraft, because the microwave power would boost the vehicle to 30000 to 50000 feet (above the 
cloud tops) - before a separate high power laser at another location was engaged. 

Fig. 49 illustrates the principle features of such a trajectory. The plot shows the pop-up 
maneuver with a moderate acceleration at liftoff which soon levels out to a constant 1 G as the vehicle 
climbs through the lower atmosphere. As the Lightcraft reaches 35000 feet (60 sec), it pitches over 
and the ERH thruster is throttled up. Near Mach 3 (90 sec) the ERH engine thrust drops rapidly 
and then the scramjet starts. Once again there is a large acceleration spike followed by a dip which 
is caused by drag build up. (The jumpy nature of the acceleration is in part due to this thrust/drag 
interaction, as well as some coarseness in the engine input data.) The sramjet “dies ” near Mach 
11 (160 sec), and the MHD fanjet is started. Acceleration during the fanjet cycle decreases as the 
coupling coefficient declines with increasing Mach number. It tails up near the end of its phase 
because the coupling coefficient levels out some (the generator exhaust is now providing a significant 
portion of the thrust), as the drag continues to decrease. Insertion is completed near 270 seconds 
and then the vehicle coasts up to the apogee (1000 sec) where a final rocket burn circularizes the 
orbit. 

The disadvantage of this latter idea is the requirement for two satellites and, of course, the 
transition between them. The high altitude satellite increases the complexity of orbital mechanics 
necessary to open a “launch window ". A traditional geostationary satellite with an orbital inclina- 
tion of zero degrees would appear quite a bit below maximum azimuth in all of the United States, 
or any location at significantly distances from the equator. The global problems of satellite power 
system architecture and relay satellite availability has yet to be studied in any detail. 

The most recent analysis is exploring the potential for combining the basic two-satellite 
boost system into a single into a single relay satellite, for fair weather launches. Fig. 50 compares 
the baseline trajectory which has a pop-up maneuver to 35000 ft. to an off-the-deck flight with a 
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Figure 50. Trajectory-Dynamic pressure. 



linearly increasing throttle setting to counter the slow pitch-over as it climbs. The dynamic pressure 
is noticeably higher in the latter case because speeds are higher at any given altitude. The total 
energy required was also about 70 GW-sec, higher due primarily to the larger drag forces on the 
vehicle. Fig. 51 shows that the primary difference in energy cost is created from about Mach 1 
to Mach 4 where the higher dynamic pressure is combined with the large drag coefficient in that 
region. After Mach 5 the curves are effectively parallel to the baseline trajectory. Another potential 
problem with off-the-deck flights is that the ERH thruster cannot reach Mach 3 until about 45000 
feet, because the ram drag is otherwise too large at these speeds and altitudes.The notch in the 
dynamic pressure (Fig. 50) is evidence that the ERH thruster is struggling to get to Mach 3. As 
the acceleration from the thruster decays the dynamic pressure actually begins to decrease before 
the scramjet kicks in. 

Much of the effort to date has been spent trying to minimize the propellant needed to 
reach final orbit and the results were surprising. Changing parameters that affect primarily the 
first part of the trajectory (i.e. pop-up altitude, ERH performance, etc.) has little effect on the 
Anal value of propellant weight. Apparently the "iterator " finds a solution so that the MHD fanjet 
portion of the flight is very similar, even though the early part of the trajectory isn’t. Changing 
performance characteristics of the fanjet or the maximum available laser power levels obviously 
affects the solution, but the general indications are that the fanjet concept can work very well in 
this type of vehicle. 

The MHD fanjet analysis section described two different generators with different assump- 
tions and performance specifications. Even when the MHD generator performance is reduced to the 
lower, more conservative generator, the results are not discouraging. The performance difference 
between the two generators gives a specific impulse reduction of a factor 1.8; e.g., if the performance 
with the optimistic generator was 9000 then the pessimistic performance is 5000 sec. This magni- 
tude of specific impulse loss would be disasterous for a chemical engine with an /, p » 400 sec and a 
propellant fraction of 85%. 

Since the propellant fraction for the Apollo lightcraft is so low (about 5-10 %), the additional 
fuel required remains fairly linear with the specific impulse change. The propellant required using 
the optimistic MHD generator is about 650 lbm and 1100 lbm for the more conservative case. This 
would keep the propellant fraction under 10 percent of the gross lift-off weight. The biggest problem 
would not be weight, but volume since liquid hydrogen has a very low density. Nevertheless, with 
such low propellant fractions the performance of this vehicle is clearly revolutionary. It is interesting 
that with propellant fractions of about 15 %, escape velocities could be attained using the MHD 
fanjet cycle. 
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The fanjet, however, is not without drawbacks. The coupling coefficient of the fanjet de- 
creases with increasing Sight Mach number, as do all air-breathing engines. At around Mach 20 
(depending on the altitude), the coupling coefficient drops below that of the rocket; therefore, pro- 
pellant weight is minimized at the expense of increased laser energy expenditure. How these two 
engine modes are to be traded against each other will depend largely on the cost of each quantity 
(i.e., LHz vs. laser energy), as well as the effects of changing vehicle configuration. 

A computer run was made to minimize the total laser energy while making very few changes 
to the input parameters.as a Fig. 51 shows the energy cost for these two cases (as well as the “off- 
the-deck " flight). The curves are almost coincident until Mach 5 where a slight difference in flight 
path angle is “lofting ” the trajectory above the baseline case - reducing drag and hence energy. 
(The switch-backs near Mach 15 are caused by an atmospheric inflection point where the pressure 
and density change abruptly; this causes the speed of sound calculation to indicate increasing values 
for a short period. There is an apparent error in the source for the atmospheric data at 170000 ft 
which causes this anomaly. The vehicle is not slowing down, but the speed of sound increase causes 
the Mach number to drop for a short period.) Near the end of ail the lines there is a horizontal line. 
This is the coast period where no energy is added. The final section shows the energy addition for 
the rocket mode. The optimized energy case clearly uses more energy in the rocket mode. Since 
the coupling coefficient is higher for the rocket than the fanjet at these speeds, the optimization 
algorithm selected a case that uses more rocket thrust than the weight optimized case. However, 
the optimized energy case uses 150 lbm more propellant than the optimized weight. 

The iterator was trying to reduce the fanjet phase (and the drag), and the only method 
available to it (in this ran) was to loft the trajectory. This causes the fanjet to terminate earlier 
because at higher flight path angles orbital mechanics dictates that the velocity be lower - to reach 
a fixed apogee altitude. This lofting is evident on Figs. 52 and 53. The primary input difference 
between the two cases is that the relay satellite is about 3 degrees higher from the horizon at pitchover 
for the lofted trajectory. Fig. 53 clearly shows the pitchover maneuver at about 60 seconds, and the 
difference in flight path angles at thrust termination. 

O. SUMMARY 

The trajectory analysis will need to continue as the engine and aerodynamic models evolve 
and the final results may end up considerably different than those presented here. This does not 
imply that the current analysis is crude or inaccurate, because it is quite precise within the confines 
of the models and inputs. The predicted performance for the Apollo Lightcraft is undoubtedly 
revolutionary, and further analysis will continue to define this new frontier. 
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Currently there are a few tender areas in the trajectory modelling effort. Among them is the 
possibility that more accurate models of the ERH and scramjet engines could leave a “hole ” in the 
thrust spectrum around Mach 3. If this was to occur another engine would be necessary (probably a 
ramjet) which would require configuration changes for the shroud and annular engine duct. Similarly, 
the MHD fanjet model does not presently have an analysis of conductivity enhancement necessary 
for significant electromagnetic power absorption. Implementation of such a model might indicate 
difficulty generating thrust at “low ” Mach numbers (i.e., near scramjet fade-out). Aerodynamic 
heating and engine heat transfer has largely been ignored and this may drive the vehicle into more 
benign flight conditions. Viscous effects (boundary layers, etc.) have been neglected and this may 
significantly alter engine performance and vehicle drag coefficients. 

Even with the above-mentioned concerns, the performance of the Apollo Lightcraft goes far 
beyond that available from any launch vehicle existing today. These early models of revolutionary 
engines are stepping stones to more powerful launch vehicles. The present analysis demonstrates 
the incredible potential of beamed-energy propulsion when united with innovative launch vehicle 
concepts. 
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CHAPTER VII 
SUMMARY 


The central objective of this NASA/USRA-supported effort was to investigate the perfor- 
mance of a revolutionary combined-cycle engine suitable for the Apollo Lightcraft. The mission 
for this advanced SSTO shuttlecraft is to transport a five person crew to low Earth orbit in three 
minutes, or anywhere on the globe in one half hour. Beam power for propulsion was limited to 10 
billion watts, a typical capacity quoted in studies of future satellite solar power stations. 

The combined-cycle engine seems to be feasible, and has now withstood a first critical 
inquiry into the propulsive physics. The group has not identified any performance projections that 
violate basic physical principles. This first-order performance analysis has unearthed a number of 
necessary refinements for each of the engine mode models, which are addressed below. 

A. ERH THRUSTER MODEL 

The operation of this innovative pulsejet engine can liken to a traditional Otto cycle com- 
plete with intake, compression, ignition (i.e combustion), and exhaust functions. A more complete 
understanding of the exhaust /intake phenomena (heretofore called “refresh ") is required to ascer- 
tain its effect upon the maximum time-averaged thrust capacity of the ERH Thruster. Also, the 
present assumptions on the peak (LSD-wave-heated) “combustion "pressure are being re-examined, 
in light of new evidence which suggest somewhat lower values. 

B. RAMJET/SCRAMJET MODEL 

Although the ramjet mode is presently “on the back burner ”, it still may be resurrected if a 
‘hole’(in propulsive thrust) appears between Mach 2.5 and 3.5. Otherwise, all attention is focused on 
scramjet engine model refinements. Since scramjet temperatures can easily exceed 10,000 R at high 
flight Mach numbers, efforts are underway to include an accurate “real gas "code. The group has a 
current copy of the NASA-Lewis “McBride "code, and plans to exercise it over the summer. Further 
model improvements are necessary in algorithms used for inlet and nozzle functions, especially with 
regard to frozen-flow effects. 


C. MHD-FANJET MODEL 
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The MHD-fanjet engine is probably the most critical “link "in the entire combined-cycle 
engine. The second most critical is the rocket mode. Unlike the ERH thruster and scramjet, they 
both have finite fuel specific impulse and dictate the total fuel load required to meet a specific Earth- 
to-Orbit mission. The ratio (by weight) of fuel consumed is roughly 80% /20% by the MHD-fanjet 
and rocket, respectively. 

As mentioned earlier, the MHD-Fanjet engine is comprised of two major components: the 
laser-heated MHD generators, and the annular MHD air accelerator. The most urgent improvement 
required in the MHD generator code is the inclusion of a credible radiative heat transfer model. 
Hydrogen gas temperatures may reach 16,000 to 20,000 K, and radiation effects cannot be ignored 
under these conditions. 

For the MHD accelerator code, a more through understanding of the air conductivity require- 
ments throughout the Mach 11-30 flight environment, is in order. Although several viable candidates 
exist for enhancing the normal air conductivity, the associated mass loss (ie., from seedant injection) 
and/or electrical power requirements must be quantified. These effects can have a significant im- 
pact on the overall specific impulse performance of the MHD-Fanjet. Furthermore, the group must 
investigate the best method for electrically uniting the generators and accelerator (ie., the power 
source and load). 

D. ROCKET MODEL 

A more detailed numerical model of the rectangular geometry LSC-wave heated Hi rocket 
gas generator is needed. This code should include non-ideal gas behavior and losses due to radiative, 
convective, and conductive heat transfer. All such effects contribute to reductions in attainable 
specific impulse. 

E. TRAJECTORY ANALYSIS 

The modified sort code is in excellent shape, but can produce results only as good as the 
engine performance models it includes. Clearly, as the engine-mode models improve, so will the 
accuracy of trajectory simulations. 

At present, the performance of the Apollo Lightcraft combined-cycle engine along an orbital 
trajectory can be described as none other than outstanding! Average vehicle acceleration is approx- 
imately 4-6 G’s, which puts the Lightcraft at orbital velocity in approximately four minutes. Beam 
power can be reduced low as 2.5 billion watts without degrading vehicle performance. Transition 


48 



between engine modes typically occurs at Mach 3, 11, and 25. The LH 2 propellant requirement for 
a 185 Km orbit is only 300Kg, or roughly 520 GW-sec for this LEO boost mission. 

If the LH 2 propellant is purchased at a price identical to that of the STS (i.e., (3.25/kg), 
the bill comes to (975; the total beamed-energy charge would total (2455, if one assumed the current 
(1.7 cents/Kw-hr that New York State pays for wholesale Canadian hydroelectric power. Therefore 
the entire round-trip-to-orbit flight aboard an Apollo Lightcraft would cost only (3430, or (686 per 
person. It is interesting to note that a round-trip charter flight to the other side of the planet by 
commercial jet aircraft (e.g. N.Y. to Tokyo) is presently (1118. Finally, it bears mentioning that 
in the Post 2020 era, space travel will become commonplace, and the present expensive “standing 
army "ground crew will be replaced by efficient macro-computers. These “tireless" workers will 
schedule launch windows, effortlessly deliver beam power from the “space power grid” to hundreds of 
spacecraft simulataneously, and electronically mail end-of-the-month billing statements to millions 
of satisfied travelers. 

In this future infrastructure, the dominant costs will be for propellant and power. Payload 
delivery costs to LEO, will finally plunge a factor of 1000 below that of the Space Shuttle Orbiter. 
Space, the Final Frontier, will no longer be a curiosity to the masses, wfc» wil mm view the Earth 
from an entirely new perspective - heretofore reserved for only a few brave astronauts in times gone 
by. 
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APPENDIX: 


THREE-DIMENSIONAL MODELING OF 


APOLLO LIGHTCRAFT LOFTLINES 
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